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Abstract 

Considered is the Schrodinger equation in a finite-dimensional space as 
an equation of mathematical physics derivable from the variational prin- 
ciple and treatable in terms of the Lagrange-Hamilton formalism. It pro- 
vides an interesting example of "mechanics" with singular Lagrangians, 
effectively treatable within the framework of Dirac formalism. We dis- 
cuss also some modified "Schrodinger" equations involving second-order 
time derivatives and introduce a kind of non-direct, non-perturbative, 
geometrically-motivated nonlinearity based on making the scalar product 
a dynamical quantity. There are some reasons to expect that this might 
be a new way of describing open dynamical systems and explaining some 
quantum "paradoxes". 

Keywords: Hamiltonian systems on manifolds of scalar products, finite-level 
quantum systems, finite-dimensional Hilbert space, Hermitian forms, scalar 
products as a dynamical variable, Schrodinger equation, Dirac formalism, es- 
sential non-perturbative nonlinearity, quantum paradoxes, conservation laws, 
GL(n,C)-invariance. 
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Introduction 

There was plenty of papers dealing with various aspects of the relationship be- 
tween classical and quantum mechanics. The most popular topics are those 
concerning the quasiclassical asymptotics of quantum mechanics, the asymp- 
totic expansions when h — > like the WKB approximation, oscillatory inte- 
grals and the method of stationary phase. There are also studies in the oppo- 
site direction, when, basing on the optical-mechanical analogy, the eikonal and 
Hamilton- Jacobi equations, one investigates purely classical structures having 
some striking analogies to quantum ones. And then one shows that really the 
mentioned classical structures, although a priori obvious on the purely classical 
level, may be also re-obtained from the quasiclassical limit transition h — > 
from quantum mechanics. Concerning such topics cf. for instance |27j and first 
of all references therein. Other very interesting studies of this kind, based on 
the hydrodynamical picture of quantum mechanics, were presented by V. V. 
Kozlov [T5]. The author analyzed there vortices of the "quantum fluid" and 
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that study is a part of his very interesting theory and methodology of vortices 
as a fundamental concept of dynamics. The ideas developed by this author go 
back to some very old and fundamental concepts in physics and philosophy of 
science. 

Another very important branch of investigations was the comparative study 
of the classical and quantum dynamics, including the dynamics of open systems 
[13l [Ml [T5] . Quantum measurement and decoherence problems in a sense belong 
to this topic [TJ [U [31 3] . Certain ideas of nonlinearity in quantum mechanics 
appeared in connection to those problems (see, e.g., [TO1 HH [Ml G2] and 
references therein). 

Our study, also as a matter of fact motivated by the aforementioned prob- 
lems, is formally quite a different approach to the quantum-classical convolution 
of concepts. Namely, we "forget" what the Schrodinger equation physically is. 
For a moment, it is for us only a differential equation. To simplify the problem 
as far as possible, we consider a "finite- level system", when the "configuration 
space" is a finite set and the corresponding linear space of "wave functions" is 
finite-dimensional. All concepts we are here dealing with will remain essentially 
valid also in the infinite-dimensional case ("true" wave functions); obviously 
some care must be taken nevertheless when passing from a finite dimension to 
the infinite one. 

Then we discuss the Lagrangian and Hamiltonian formalisms for such a 
"classical" mechanical system in C". In particular, the Dirac theory of primary 
and secondary constraints for systems with degenerate Lagrangians is discussed. 
This approach enables one to formulate some models of nonlinearity. We hope 
this nonlinearity may be perhaps a tool for describing measurement paradoxes 
and decoherence. The main point is that our formalism seems to suggest in a 
natural way some geometrically well-justified nonlinearities, not ones introduced 
"by hand" as perturbations of some linear background. The main idea is that of 
dynamical " scalar product" which is not fixed once for all but itself is a dynam- 
ical quantity on equal footing with the "wave functions"; they both satisfy a 
closed system of essentially, non-pcrturbatively nonlinear differential equations. 
The structure of this nonlinearity is based and, one can say, almost canonically 
implied by the geometric structure of "classical" degrees of freedom. Because 
of this we hope that the resulting effectively nonlinear quantum mechanics may 
be perhaps free of paradoxes of decoherence and measurement and can provide 
some new description of open quantum systems, alternative to that described 
in [DEI El . 

Some ideas of our "classical" description of quantum systems in terms of 
phase spaces and Hamiltonian dynamics are similar to those suggested many 
years ago by D. Chrusciriski [S]. 
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1 Finite-level nonlinear Schrodinger equation in 
the Lagrange-Hamilton description 

Following the jargon used by laser specialists and those working with the quan- 
tum dynamics of mutually interacting spins, we use the term "finite- level quan- 
tum system" for such a one the " Hilbert space" of which is finite-dimensional, so 
it may be identified with C", when some basis is fixed. However, we shall avoid 
the misuse of this identification, because it usually smuggles into the treatment 
some artificial objects obscuring and often just falsifying the proper geometric 
interpretation of the used concept and making impossible the introduction of 
new ideas. 

And one thing must be explicitly stressed. We are in fact motivated by 
certain problems from the realm of foundations of quanta. However, for some 
reasons at this stage of our treatment it was convenient to pretend to "forget" 
about this motivation and just to consider Schrodinger equation in a finite- 
dimensional space purely classically, simply as an equation of mathematical 
physics, in a sense classical mechanics, derivable from the variational principle 
and because of this treatable in terms of Lagrange-Hamilton formalism. And it 
is really interesting even from the point of view of this "hypocritically" classi- 
cal language. For instance, it provides an interesting example of "mechanics" 
with singular Lagrangians, effectively treatable within the framework of Dirac 
formalism involving the primary and secondary constraints in a phase space of 
the problem. 

Later on we shall try to discuss the " Schrodinger" equation involving second- 
order time derivatives, and also introduce some kind of non-direct and geo- 
metrically-motivated nonlinearity based on making the scalar product a dy- 
namical quantity. There are some reasons to expect this might be a new way 
of describing open quantum systems and a new promising attempt towards ex- 
plaining quantum "paradoxes", decoherence and measurement. And using the 
methods developed for n-level quantum systems, we formulate finally some ideas 
concerning the treatment in an infinite-dimensional Hilbert space and certain 
links with relativistic field equations. 

1.1 Some complex geometry in linear spaces 

In this section our primary concept is an n(< oo)-dimensional linear space W 
over the complex field C. It is well known that such a space gives rise to the 
natural quadruple of mutually related n-dimensional complex linear spaces: 

w, w*, W, W* ~w^, (i) 

namely, W itself, its usual dual W* over C, i.e., the space of C-linear C- valued 
functions on W, and their complex-conjugate spaces W and W ~W*. There 
are many mistakes and misunderstandings concerning the complex conjugate 
space W and the antidual one W ~ W* 7 which may be easily avoided at 
least if W is finite-dimensional, that is the case at the present stage. Hence, 
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W ~ W* consists by definition of antilinear (semi-linear) functions on W, i.e., 
such ones which satisfy 

g(aw + bv) — ag(w) + bg(v) (2) 

for arbitrary w, v € W , a, b S C. So, there is a natural antilinear isomorphism 
of W* onto W* given by 



W~ 3 /■-►/€ W*, f(w):=f(w) (3) 

for any id £ W, f £ W*; the complex conjugate is taken pointwisely. As 
mentioned, ([3]) is an antilinear isomorphism acting between two linear spaces, 



(af + bg) = af + bg (4) 

for any /, g € W* , a,b E C. For obvious reasons the inverse of ([3]) will be 
denoted by the same symbol and the following holds: 

7=/- (5) 

By analogy to the obvious canonical isomorphism between W and W** , 
the complex-conjugate space W is defined as the antidual of W* . Namely, its 
elements u € W are by definition antilinear functions on W*. So, for any 
feW*, u(f) is defined as 



u(f) := f(u) = f(u). (6) 

Compare this with the standard identification ofw£W with the linear function 
on W*: 

w(f) = f(w). (7) 
Again there exists a canonical antilinear isomorphism of W onto W: 



W 3 w ^ w S W, w(f)=w(f). (8) 

All this may be also done in other, equivalent way. The important things are 
the following ones. All linear spaces (Q} are logically different and the complex 
conjugates of vectors belong to other linear spaces, for example, W £ W is not 
an element of W , unless some additional geometric structures are fixed in W. 
There exists canonical linear isomorphism between W and W** (and similarly 
between W and W**). There exists also canonical antilinear isomorphism of W 
and W (and similarly between W* and W*). But of course without additional 
geometric object in W (some metric) there is neither any canonical isomorphism 
between W and W* nor any one between W and W* . 

If (ei, . . . , e„) and (e 1 , . . . , e") are mutually dual bases in W and W* , 

e a (e b ) = 6\, (9) 

then by the complex conjugation one obtains from them canonically the mutu- 
ally dual bases (ej,... ,e«r), (e 1 , . . . ,e") respectively in W and W ~ W*, 

^0%) = 5\ (10) 
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If we expand vectors and covectors in W with respect to the first system of dual 
bases, 

w = w a e a , f = f a e a , f(w) = f a w a , (11) 

then, obviously, 

w = vfea, 1 = l^t. (12) 

In this special system of bases the complex conjugation of objects is analytically 
represented by the complex conjugate of the coefficients, as seen from the above 
formulae. But certainly w, w, /, / are elements of pairwise different linear 
spaces. So, the above formulae introduce the system of notations which will be 
consequently used below. 

Let us stress that the definition of W* given above is valid in any linear space 
independently of its dimension being finite or not. Unlike this the definition of 
W quoted above works only in the algebraically reflexive linear ones, but no 
other ones present interest for us. 

Further on, the whole system of tensor products of W, W* , W, W* may be 
introduced. Obviously, in applications we are interested in, the most important 
are linear mappings of W into itself, i.e., elements of 

L(W) ~W ®W* (13) 

and sesquilinear forms, i.e., elements of W*®W* . The corresponding analytical 
matrix representation is given respectively by L a b, F^b for L e L(VF), F e 
W* <g> W*, etc. For sesquilinear forms we use the convention that they are 
antilinear in the first argument and linear in the second one: 

F(u,v) = FabtFv b , (14) 

so that 

F(au+bw, v) — aF(u, v)+bF(w, v), F(v,au+bw) = aF(v,u)+bF(v,w) (15) 

for any values of vectors and coefficients. We usually need Hermitian forms, i.e., 

such ones that 

F(u, v) = F(v,u), Fab = F b a. (16) 

Then, obviously, 

F(u,u) = F(u,u), u^F(u,u)eM.. (17) 
If F is non-degenerate, 

det [Fa b ] ^ 0, (18) 

i.e., if 

W 3 v ^ F(-,v) e W 7 (19) 

is a linear isomorphism of W onto W* 7 then there exists the inverse twice 
contravariant sesquilinear object F^ 1 e VF®VF with components denoted briefly 
by F al such that 

F a ~ c F- cb - S\, F Uc F c ~ b = 5a». (20) 
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Obviously, S a b, 5a are here matrices of identity transformations respectively in 
W and W 7 . 

If F is Hermitian and for any u e W we have 

F{u,u)>0 if u^O, (21) 

then we say that F is positively definite. Then there exist such bases in W that 
the matrix of F is diagonal with (+l)-entries on the diagonal, 

[Fab] = In- (22) 

Otherwise F has some signature (the number of positive and negative elements 
on the diagonal) which is an invariant of F. 

Let us now assume that some non-degenerate Hermitian form T G W* ® W* 
is fixed once for all in W, so we are dealing with the algebraic structure (W, T). If 
r is positively definite, (W, T) is a finite-dimensional Hilbcrt space or, using more 
customary terms, a unitary space. If the positive definiteness is not assumed, 
we are dealing with a pseudounitary space. The form T will play a role of 
the scalar product if we insist on the quantum-mechanical interpretation. So, 
it seems natural to expect it should be positively definite. However, in many 
formal problems this assumption is not necessary. Moreover, even within the 
framework of quantum-mechanical interpretation, one cannot exclude a priori 
models with non-definite scalar products, at least at some stage of considerations 
(remind the old quantization of electrodynamics according to the Gupta-Bleuler 
prescription) . 



1.2 Lagrangian and Hamiltonian models of the usual and 
generalized Schrodinger equations (evolution) 

Let us return to our " classical hypocrisy" . We discuss some models of " analytical 
mechanics" in the "configuration space" W. The system has 2n degrees of 
freedom because the dimension of If as a linear space over reals equals 2n. 
However, often it is convenient to say that we are dealing with n "complex 
degrees of freedom". If some basis (ei, . . . , e n ) is fixed in W, then expanding 
the elements \p <G W, 

V> = 4> a e a , (23) 

we introduce "complex generalized coordinates" i[) a . Their real and imaginary 
parts form a system of 2n usual real coordinates. It is more convenient to use 
another normalisation, namely, 

r = ^(x a +*y a ), ^ =^={x a -iy a ), (24) 

where x a ,y a e M. Obviously, without any additional structure in W, the real 
and imaginary parts of the vector ip <G W are not well defined. Taking the real 
and imaginary parts of ip a is an explicitly base-dependent procedure. Just like 
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in certain well-known formulae of classical field theory it is often formally conve- 
nient to use the system of 2n complex variable ip a , ip as "coordinates". Many 
expressions and calculations look then simpler. Then 2n complex quantities tj) a , 
iff 1 are (on some intermediate stages) treated as if they were formally indepen- 
dent and only in the final formulae one " remembers" that they are interrelated. 
And expressions (|24[) describe then something like an orthogonal transforma- 
tion of the system of 2n real variables (x a ,y a ) into the system of 2n complex 

variables (^ip a More precisely, in a sense we work in the complex space 

W x W of the complex dimension 2n. And later on the formulae are restricted 
to the "diagonal" of complex dimension n (real dimension 2n) consisting not of 
all pairs (ip, ip) € W x W but only of the pairs of the form (ip, tfjj <E W x W . If 
there is some need of using symbols, we shall denote this "diagonal" as 

Diag^xW 7 ) :={(^,V0 -ipeW}. (25) 

In variational problems one deals with real- valued Lagrangians. It is convenient 
to define them primarily as analytic functions on W x W and perform all dif- 
ferential operations with respect to both arguments as independent ones. Of 
course, if non-constant, such functions are never real-valued. However, they 
are constructed so that to be real- valued on Diag (W x W) . If they are ex- 
panded into double power series with respect to the (W, W) -arguments, their 
coefficients must show Hermitian symmetry. To be more precise, this concerns 
functions of the potential energy type, i.e., which depend only on configurations. 
Lagrangians, however, depend on configurations (generalized coordinates) and 
velocities. So, they are functions on WxWxWxW~WxWxWx W. 
Nevertheless the above statements concerning structure of potentials apply also 
to such tangent-bundle functions: simply W, W are then replaced by W x W, 
W xW ~ W x W. 

1.2.1 Lagrangian for the standard Schrodinger equation 

Let us begin with the usual Schrodinger Lagrangian. Later on we shall consider 
the hierarchy of more and more complicated Lagrangians with the additional 
terms responsible for the various expected physical phenomena. 
Definition: The usual Schrodinger Lagrangian will be denoted by L(l), where 
the label (1) refers to the resulting first-order differential Euler-Lagrange equa- 
tion. It is given by 

L(l) = iar m (VV - ^V) - lXabTi>\ (26) 

where a, 7 are some constants, \ E W* <g> W* is a Hermitian form on W, and 
obviously the dot-symbols are time derivatives. 

It is interesting to mention that Lagrangians of this form appeared in a 
slightly different context in [23l [26] . They were thought on as an alternative 
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variational description of vibrations. And besides, the exchange of energy be- 
tween real and imaginary parts of i\) a was expected to be some toy model of 
dissipative phenomena (the observed system and its unobserved partner) . This 
has to do with so-called Birchhoffian formulation of dynamical laws. 

Similar expressions appear also in Kozlov theory of vortices |19j and in his 
study of optical, quantum-mechanical and quasiclassical phenomena. In partic- 
ular, this concerns vortices appearing in hydrodynamical interpretation of the 
Schrodinger equation and its quasiclassical limit. 

Due to the hermiticity of T and %, L(l) is real, i.e., 



L(1) = L(1), (27) 
just as Lagrangian should be. The corresponding action is given by 

1(1) = { L(l)dt (28) 



2iaTa b ip b - lXab^\ (29) 



and its variational derivative equals 

(5/(1) _ dL(l) d dL(l) 



where, according to the standard procedure mentioned above, ip and ip are 
formally treated as independent quantities. Then, obviously, calculated in the 
dual way 

5ij} a (t) dip a dt dip a 

is the complex conjugate of (|2"9"|) . Hence, it is only one of these expressions, 
by convention (|29[) . that is sufficient for obtaining equations of motion (Euler- 
Lagrangc equations). Putting (|29p to vanish one obtains after some manipula- 
tions on tensor indices the following equation: 



>>a— = - ( X) bi> , (31) 



where T x — H£ Lc(VK) is obtained from x G W* ® W* by the T-raising of the 
first tensor index, i.e., 

H\ = ( r X )\ = T' rc x- cb . (32) 

Obviously, H — T x is by its very structure Hermitian with respect to T (x itself 
is objectively Hermitian), 

T(H1>,v) = T{1>,H<p). (33) 

It is seen that (|3"Tj) becomes the literally understood Schrodinger equation 
when r is positively definite, so in appropriate bases 

[Tab] = In (34) 
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and a = h, 7 = 2, where, obviously, h is the Planck constant. 

Lagrangian (|26p is linear in generalized velocities ip, ip with coefficients de- 
pending linearly on generalized coordinates ip, ip. Therefore, the action func- 
tional (f28|) is quadratic with respect to the evolution curve 1 9 t h> ^(t) 6 M / 
and its " stationarization" results in first-order linear differential equations for 
the time-dependence of ip. As a matter of fact x? an d then also if, may be 
time-dependent, and then the solution is not simply given by the correspond- 
ing operator exponent. In "usual" analytical mechanics Lagrangians linear in 
velocities and first-order differential equations of motion are rather exotic (al- 
though second-order differential equations may be in an obvious way reduced 
to doubled systems of first-order equations). And the resulting Legendre trans- 
formation is evidently non-invertible, one obtains constraints in the phase space 
of the system and the Dirac formalism of singular Hamiltonian mechanics must 
be used. As we shall see, this formalism for the Schrodinger equation has some 
geometric peculiarities and is interesting in itself. 

1.2.2 Admitting second derivatives 

If we continue to forget "hypocritically" about our quantum motivation, then 
from the point of view of purely classical analytical mechanics it is natural to 
ask after Lagrangians which lead to terms with second derivatives in equations 
of motion. 

Definition: Let us denote those Lagrangians by L(l, 2), where the labels (1, 2) 
refer to the occurrence of first and second time derivatives (velocities and accel- 
erations) in equations of motion. Obviously, the simplest model is 

L(l, 2) = laTab (VV b - + /OTWV 6 - IXabT^- (35) 

If a = 0, this becomes just the usual Lagrangian of the system of 2n coupled 
harmonic oscillators. However, being motivated by the quantum-mechanical 
problems and the usual Schrodinger equation, we are inclined to retain the 
a-term; then the quadratic /3-term is some kind of correction. Its physical 
interpretation in quantum-mechanical terms is still not clear, if possible at all. 
Remark: One might formally admit to take another T in the term quadratic 
in velocities, then Lagrangian would have the following form: 

L(l, 2) = L (l, L; 2, f ) = iaT m fa?ip b - + /9IW > '-Ttttt? V '• (36) 

Of course, this makes the resulting equations much more complicated and one 
feels rather reluctant to such a modification, nevertheless it is formally possible. 
And perhaps it may be physically justified, provided of course that the quadratic 
correction may be physically interpretable at all. 

Let us observe that from the point of view of purely classical analytical 
mechanics the quantities P, L are logically independent and it is fully justified 
to discuss dynamical models in which they are different and non-correlated to 
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each other. Obviously, the resulting equations of motion would be then more 
complicated and the symmetry group rather restricted (because it must preserve 
two different sesquilinear forms P, T). The mathematical peculiarity of P, T 
proportional to each other is just the "large" symmetry group of equations of 
motion. 

Nevertheless, our analysis has some "quantum" motivation, where T is to 
play the role of scalar product (or some generalised scalar product). And then 
it is more reasonable and physically justified to use only one form T. 
Remark: Structural similarity of Lagrangian (|26p to that underlying Dirac 
equation is obvious, although Dirac equation is a partial one, and here for sim- 
plicity we deal with ordinary differential equations (finite " configuration space" ) . 
This similarity is not accidental (see e.g. our papers [551 UHl GUI OS] ) ■ The same 
concerns (|35[) where in the resulting differential equations second-order time 
derivatives are combined with first-order ones [20l [22] . 

The corresponding action functionals will be denoted respectively by 

1(1,2) = f L(l,2)dt, 7(1,2)= Jl(l,2)dt. (37) 

Putting a = we obtain the "usual" Lagrangian of analytical mechanics (for 
2n coupled harmonic oscillators, as said above). The corresponding Lagrangian 
will be denoted by L(2), 

L(2) = /HwV ~ lXmTi>\ (38) 
and its action functional by 1(2), 

1(2) = J L(2)dt. (39) 

It is obvious that 

— = — = 2iaT m> — - P^ab-j-r- ~ JXab4> ■ (40) 
Sip (t) dt dt 

There is a qualitative catastrophic discontinuity if performing the limit transi- 
tion P — > in the corresponding equations of motion: 

ia^ - = lH^, (41) 

dt 2 dt 2 2 ^ ' V 1 

where the "quantum" Hamiltonian H a b is given again by (|32p . This is nothing 
else but the special case of the general phenomenon that differential equations 
and phase portraits of the corresponding dynamical systems are drastically un- 
stable with respect to neglecting the highest-order derivative terms. 

Obviously, on the level of pure analytical mechanics it is just (|38p and ([3"9"]) , 
i.e., a = situation, that is the most natural model. If we do not hide the (true) 
quantum-mechanical motivation, then obviously the model (|26p . i.e., j3 = sit- 
uation, seems to be just the true model. The natural question arises as to the 
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status of the second derivatives. Certain arguments were raised by various au- 
thors in favour of the /3-term in L(l, 2), i.e., "quantum mechanics" corrected by 
the second- derivative ("acceleration") expression in the equation of evolution 
for the "wave function" tp [TTJ [T71 [2TJ 131] ■ It does not seem yet clear, how- 
ever, if this term may be made compatible with the statistical interpretation. 
And if it is incompatible, what might be a possible alternative interpretation to 
be used instead. This finite-level problem resembles the interplay of first- and 
second-order derivatives of wave amplitudes in Dirac and Klein-Gordon relativis- 
tic equations. And it is really a " discretized" version of the problem. Moreover, 
the mixing of terms like in (f3"5)) and (|4"Tj) is known from the field-theoretic treat- 
ments. There exist relativistic models in which Dirac and d'Alembert operators 
are superposed [Til HZ1 HU [HI [29l [301 EU • There were various motivations for 
that. One of them was the demand of conformal invariance and the gauge model 
of gravitation based on the conformal group [17l [28l [29j [301 E2] • Incidentally, 
the mentioned models have some quite unexpected and interesting consequences 
which may be of some relevance for the standard model. 

We do not discuss here the model (f3T>)) with different Hermitian metrics V, T 
for the first- and second-order terms. Nevertheless, let us quote the expression 
for its variational derivative and the corresponding "Schrodinger equation" in 
the form solved with respect to the first derivatives: 

dib b ~ d 2 tb b 

2iaTm ^T ~ pTm ^k ~ 1Xm%pb ' (42) 

/?( r r) b ^ +1 H\^\ (43) 
where T T £ L(W) is given by 

( r f)\ :^T a ~f- cb (44) 
and H = T \ G L(W) ~ W ® W* is, as usual, the " Hamiltonian" (gg) . 

1.3 Admitting "potentials" and direct nonlinearity 

The above models are linear and finite-dimensional with perhaps some time- 
dependent coefficients (x, or equivalently H). Therefore, they are formally 
equivalent to a finite system of harmonic oscillators with a possible parametric- 
type excitation. Incidentally, to include non-parametric-type excitation, one 
has to admit Lagrangians to be general second-order polynomials of the state 
variables, not necessarily quadratic forms. The terms like 

icxc = F a 4> a + FaT (45) 

with time-dependent F a describe in analytical mechanics the extra-imposed ex- 
ternal extortion. The resulting Euler-Lagrange equations are affine (" linear- non- 
homogeneous" ) but no longer literally linear in state variables. So, although nat- 
ural and very well known in analytical mechanics, they are outside the scope of 



dlp a 
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quantum mechanics with its linear Schrodinger equation. But being once faced 
with such an "elementary" nonlinearity ("pseudo-nonlinearity" , so to speak, in 
the sense " linear- non-homogeneous" or "afhne"), one feels motivated to admit 
a real, serious nonlinearity by introducing to L some general potential. In ana- 
lytical mechanics this is something very natural and belonging to the everyday 
practice. The "usual" quantum mechanics is linear and within this framework 
such corrections might seem exotic. Nevertheless, it is well known that there 
exists some motivation and there were some attempts to introduce nonlinearity 
to quantum mechanics. They have to do with so-called "paradoxes" like the 
reduction of wave packet, decoherence, measurement, etc. Many of such at- 
tempts were "blind" in the sense of introducing nonlinearity "by hand" using 
the trial-and-error methodology. The language of geometric models in analyti- 
cal mechanics enables one to proceed in a more systematic way basing on some 
natural guiding hints. 

As mentioned, the simplest way is to introduce to Lagrangian some potential 
term V built in a non-quadratic way of the "wave function" ip. This is still a 
" traditional" way based on some kind of ad hoc " experiments" with postulated 
potentials. The method we suggest in this paper, namely the one based on the 
dynamical scalar product, seems much more natural and promising. Neverthe- 
less, let us mention in a few words the traditional procedure. 
Definition: The corresponding Lagrangians will be denoted by L(dnl), where 
"dnl" means "directly nonlinear" . They may be based cither on the Schrodinger 
linear background L(l) (f26|) or on the linear background L(l, 2) (|35|) predicting 
the second- derivative term, 

L(l,dnl) = L(l) + V, (46) 
i(l,2,dnl) = L(l,2) + V, (47) 

where V is just the potential term built in a non-quadratic way of ip and re- 
sponsible for the " direct nonlinearity" . 

The corresponding contribution to the action functional will be denoted by 

I(V) := J Vdt. (48) 

If we use the complex formalism, it is convenient to follow the procedure de- 
scribed above. Namely, we start from some analytic function of two variables, 
i.e., V(ip,<p), ip € W, ip G W, satisfying the aforementioned hcrmiticity condi- 
tion, therefore real on the diagonal {(ip,^) : ip € W}. 

Proposition: The most natural expressions are ones built of the obvious in- 
variant Tabfp "0 b j 

^ (VS^ = / (WV) , (49) 

where / : R — > M is some real-valued function of one real variable. 

For example, one can think about the "quartic" model often used in quantum 
field theory and elementary particles physics: 

f(x) = k(x - a) 2 , (50) 
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where k, a £ R denote some real constants. Obviously, the variational derivative 
of the corresponding functional I(V) is given by 



5 Sp- = f (iWV) (5i) 

Sip (t) v 1 

where /' is the usual first-order derivative of / : R — ► R. 

Variation with respect to ip a is given by the complex-conjugate expression: 



8ip a {t) 



!' (iWV) ip b T- ba . (52) 



The corresponding nonlinear Schrodinger equation with the possible second- 
order differential term has the following form: 



dib b d 2 ib b 
2iaT m ^- - l3Ta b -^r = lXm4> b + /T a f, (53) 
at dt z 



i.e., 



l a d -^-^ = lH^ + l -rV. (54) 
dt 2 dt 2 2 v 2 J v v ' 

1.4 Canonical formalism 

We shall now discuss some problems of canonical formalism for the above usual 
and modified Schrodinger equations. Before doing this we again return to some 
comments concerning our complex language. As mentioned, just like in some 
studies concerning classical field theory and its quantization, it is convenient to 
use the complex formalism. Let us adapt it to the phase-space description. As 
mentioned, in a fixed basis in W we put 

r = {* a + W a ) , T = -^=(x a - iy a ) , (55) 

where x a , y a form a system of 2n real coordinates in the configuration space. 
In other words, we perform analytical continuation from the real form of W 
to the complex space W x W and then perform the restriction to the diagonal 
Diag(W / x W) (USD. If we use the language of real geometry, the canonical 
momenta conjugate to (x a ,y a ) are denoted respectively by (u a ,v a )- They are 
coordinates in the dual space W* (as a real space). In the language of complex 
geometry, we use the momentum space W* x W , i.e., the phase space W x 
W x W* x W*. The momenta conjugate to (^ip a ,ip a ^ are denoted by (7r a ,7Fa), 
where obviously 

7T a = — i= (U a - iv a ) , 7Tq = -j= (u a + W a ) . (56) 
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Inverting formulae (|55|) . (|56")l we obtain 
1 / 



I 

71 



(V - / 



71 ( " a " ' 



(57) 
(58) 



Formally we are dealing here with the orthogonal change of variables from (x, y) 
to (ip,^p) with the block matrix 



1 

V2 



Hn 

-il„ 



(59) 



where /„ obviously denotes the n x n identity matrix. Transformation from 
(u, v) to (it, 7f) is obviously given by the contragradient (inverse and transposed) 
matrix. The relationship between differential operators in (x,y;u,v) and the 
corresponding analytical continuation to (ip, tp; tt, 7t) is evidently given by 



1 d 



1 d 



a 

did Id 

du a y^dlTa ^dWa' 

and conversely 



d 



d 



d 

dy a ~ ^2diP a V2a^' 
d id id 

dv a V2dn a ^/2dTfa' 



d 


1 d 


i d 


d 


1 d 


i d 


dip a 


y/2dx a 


~ V2dy°-' 


w 


^/2dx a 




d 


1 d 


i d 


d 


1 d 


i d 


dlTa 


~ V2du a 


V2dv a ' 


9Wa 


~ V2du a 


V2dv a 



The symplectic form 

7 = du a A dx a + dv a A dy a , 
after analytical continuation looks as follows: 

7 = dn a A di/j a + dWa A dip a , 
and, as expected, the Poisson brackets 

df dg df dg df dg df dg 



{f,g} = 



dx a du a dy a dv a du a dx a dv a dy a 
after analytical continuation become 

df dg df dg df dg df dg 



if, 9} 



dtp a dn a Q^ a dWa dn a dtp a dWa Q^> 

Hamiltonian vector fields with generators F, 

dF d dF d dF d dF d 

J£ p, _j_ 

du a dx a dv a dy a dx a du a dy a dv a ' 



(60) 
(61) 

(62) 
(63) 

(64) 
(65) 

(66) 
(67) 

(68) 
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become obviously 

x _ dF d dF d dF d OF d 
F dn a dip a dW-aQ^ dip a d-K a g^jf dWa' 

Let us quote some additional obvious formulae which often appear in our 
calculus concerning the above and other formulae: 




Similarly, 



d^)=S a \ (^»^)=0, (73) 

Concerning the formulae like (|65|) . (|67|) . (|69]) . and so on, it must be stated 
that working in the "configuration space" W x W and the "phase space" W x 
W x W* x W is an auxiliary tool, although very convenient one. However, 
the true "physical" phenomena take place in W and W x W* as the configu- 
ration and phase spaces, respectively. Fortunately, the analytical continuation 
from Diag (W x W) enables one to work almost automatically in the mentioned 
spaces of the doubled dimension. 



1.4.1 Legendre transformation, constraints, Dirac procedure 

Let us consider the Legendre transformation based on the usual Schrodinger 
Lagrangian (f!?6")) . As the corresponding action functional if) i— > is quad- 

ratic and the resulting "equations of motion" are linear, this study is rather 
academic. Nevertheless, it remains essentially valid when some direct nonlin- 
earity is introduced by the "potential energy" term V (if), ip) and everything 
becomes especially instructive when the term quadratic in velocities is intro- 
duced, i.e., if we consider the model £(1,2) (f3"5| and the problem of the limit 
transition f3 ~ * 0, a very singular one. For the pure model L(l) the Legendre 
transformation has the following form: 

n a = iaT- Sa ip b , Wa = -iaTabip h ', (74) 

the second equation being the complex conjugate of the first one; this is the 
obvious redundancy following from the use of complex language. It is seen 
that canonical momenta are completely independent of generalized velocities, 
so this is an extreme case of singular Dirac mechanics (compare the situation, 
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by the way, with one for the relativistic Dirac equation). Therefore, the primary 
constraints are redundantly described by equations 



a =0, <fe = 0, (75) 

where 

<f>a = TTa ~ io^^tp , (j)^ = 7T „ + iaTabi> h '. (76) 

According to the general rules of Lagrangian-Hamiltonian mechanics, the "en- 
ergy" function e is given by 

in this case L = L(l) ([26| . Of course, one should not confuse this "energy" in 
the sense of Hamiltonian mechanics with quantum energy given by the Hamilton 
operator (f52| . After some trivial calculations we obtain that 

e - 7Xafe?V, (78) 

so the "energy" function e is independent of generalized velocities — a rather 
exotic property. As usually, e is a pull-back, under Legendre transformation C, 
of some " Hamiltonian" h defined only on the manifold 

M = C (W x W x W* x W*) (79) 

of primary Dirac constraints in the phase space of the system. However, accord- 
ing to the traditional procedure of Dirac, it is convenient to use the family of 
" Hamiltonians" H defined on the "total phase space" W x W x W* x W and 
having the property that 

H\M = h. (80) 
The standard procedure is to fix any such a Hamiltonian H.q and then to put 

H=H Q + \ a ^a+T^a (81) 

with yet undetermined Lagrange factors A a , A . The most natural, almost 
canonical choice is _ 

Ho = ixmT^- (82) 

Then, following the well-known Dirac procedure, one must determine the sub- 
manifold M s of secondary constraints, eliminate as fast as possible the above 

Lagrange coefficients A a , A a ("gauge variables"), define the effective Hamilto- 
nian on M s and the corresponding Poisson brackets of functions on M s (Dirac 
brackets in a sense). The first step is to calculate Hamiltonian vector fields X^, 

Xh = X Ha + X\a,p a + Xj--^_. (83) 
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More precisely, it is sufficient to take the simplified form of Xn, namely, 

X H = X Ho + \ a X^ a + 1\, (84) 

because the expressions ([83]). t|84|) evidently coincide on the manifold of primary 
constraints. Then, the care must be taken to make X-h compatible with con- 
straints M , i.e., being tangent to it. At points of non-tangency the dynamics is 
inconsistent. 

It is a trivial task to calculate the contractions of differentials of (|75|) . (I76| 
with the vector fields |83|) . (|84|) . As differentials of constraints equations are 
given in the redundant space W x W x W* x W by 

d(f> c = dn c — iaT- lc dip d , dcj)^ — dn- + iaTcddi{> d , (85) 
the tangency conditions 

(# C) X H ) = 0, (<% X H ) = (86) 
have obviously the unique solutions for A a , A , i.e., 

A a = -^r»W b , T = ^^ c r ffi , (87) 

all over the primary constraints manifold M. No additional restrictions are 
imposed on the admissible points of M and because of this the secondary con- 
straints M s (in the Dirac sense) coincide with the manifold of primary con- 
straints. No gauge freedom appears because at all points of M = M s the La- 
grange multipliers are uniquely defined by ([87)1 . The momentum variables ir a , if a 

arc uniquely determined by the generalized coordinates ip a , ip . Let %m = im„ 
denote the natural injection of M — M s into the phase space manifold, then 

7 ||M = 7 ||M S =^7 = ^7, (88) 

i.e., the restriction of the natural symplectic two- form on the primary phase 
space P to M = M s is simply given by 

7 | \M = 7 | \M S = 2iaT 1Ib <0 : A # b . (89) 

Roughly speaking, this means that on M = M s the quantities -0° become ef- 
fectively the canonical momenta conjugate to tf> a as the generalized coordinates 
(obviously, up to the linear transformation with matrix 2iaTab)- 

Some subtle points appear here due to the use of complex coordinates ip a - 
Since their number is arbitrary, it takes any admissible value n, not necessarily 
an even one as it must be in symplectic manifolds. But M — M s has the 
complex dimension n, whereas its real one is always even and equals 2n, cf. the 
formulae ([55)1 - ([55]) . (f6T)]) - (|6"5"l) . Therefore, the n complex parameters 

II a = 2ia^ b T- ba (90) 
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really provide the complex representation of the canonical momenta conjugate 
to tp a . Let us describe this effectively in terms of real parameters (f57|) . (|58|). 

Some remark is necessary here to avoid the conflict (the multiplier 2) between 
the last comment concerning and the formulae (f?4")) describing the Legendre 
transformation in terms of complex variables. Namely, (|59"| and (|9T)|) describe 
the effective Darboux representation in the constraints manifold M = M s , not 
in the original non-restricted phase space. 

Take the formula (|89[) and express it in terms of some fixed basis (. . . , e a , . . .) 
introducing the corresponding real coordinates x a , y a (|5"T|) . Let us express Tab 
in terms of this particular basis, 

Tah = S ab + iA ab , (91) 

where S ab , A ab are respectively symmetric and anti-symmetric real matrices, 

Sab = Sba, A ab — — A ba , (92) 

always related to this particular choice of basis. Then after straightforward 
calculation one obtains 

7 ||M = 7 ||M S = -2aS ab dx a A dy b - uA ab (dx a A dx b + dy a A dy b ) . (93) 

Proposition: If the complex basis (. . . , e a , ■ ■ ■) is chosen in such a way that 

A ab = 0, Sab = ^-S ab (94) 

(it is always possible when T is Hermitian and positively definite), then 

7 ||M = 7||M S = S ab dy b A dx a = dy a A dx a . (95) 

So, we conclude that indeed, from the point of view of real linear structure, 
M = M s is a real symplectic manifold with Darboux coordinates x a (effective 
generalized coordinates) and y a — 5 a by b (generalized conjugate momenta). So 
indeed, from the point of view of real symplectic geometry, M — M s is a purely 
second-class manifold (in Dirac language). The effective Hamiltonian responsi- 
ble for the L-dynamics is given on M = M s by 

H = H \M = H \M S , (96) 

i.e., ([52")) . Let us substitute again the coordinates as above and put 

Xab = Oab + ia a b, (97) 

where again c, a are real and respectively symmetric and anti-symmetric ma- 
trices. We obtain 

H = \crab (y a y b + x a x b ) + ^a ab (x b y a - x a y b ) . (98) 
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This is the real form of the reduced Hamiltonian of dynamics derived from the 
Schrodinger Lagrangian. The equations of motion have the form of the Hamilton 
equations with the symplectic form (|95p . x a , y a — 6 a bU b being respectively 
canonical coordinates and their conjugate momenta, with the standard Poisson 
brackets. 

Remark: The first term in (|9"5)) is the usual isotropic harmonic oscillator 
(isotropic in the sense of a a b, not S a t, unless both are proportional to each 
other); the second one is not particularly important because for the Hermitian 
bilinear forms Xab, ^sb, when in addition T is positively definite, it is always 
possible to choose a basis (. . . , e a , . . .) in which simultaneously T is 5-like and \ 
is symmetric, i.e., 

a ab = 0, (99) 

moreover, x is then real-diagonal. Formally it is not necessary to assume that 
[Xab], [&ab] must be positively definite. If they are not, the mentioned interpre- 
tation in terms of the usual harmonic oscillator is not literally true. 

It is perhaps instructive to express everything in terms of the real variables 
(|57|) , (|58|) . although the complex ones (|55|) , (|56|) are formally more convenient, 
even in field-theoretic problems. We did this partially in formulae (|64|) — (|69|) . etc. 
Let us now look what is the real representation of the Legendre transformation 
(|?4")) and some resulting relationships. It is easy to show that the Legendre 
transformations become 

u a = aA ab x b + aS a bV b , v a = -aS ab x b + aA ab y b . (100) 

Proposition: If the basis (. . . , e Q , . . .) is chosen in such a way that (|94|) holds, 
then one obtains simply that 

U a = ^S ab y b , V a = -^6abX b . (101) 

Proposition: It is perhaps a little pretentious and artificially sophisticated, but 
nevertheless instructive and aesthetic to admit more general bases (. . . , e a , . . .) 
in which 

Sab = ^gab, (102) 

la 

where g a b = g ba and in "physical" models the matrix [g a b] is positively definite. 
It plays the role of the Euclidean metric tensor in the n-dimensional real linear 
space U composed of linear combinations X a e a with real coefficients A a e K 
(the R-linear shell of the system of vectors e a , a — 1, n). Obviously, the linear 
quantities 

Va ~ g ab y b (103) 

may be interpreted as components of t/-covariant vectors, y € U*. Then (|93p 
becomes 

7 | \M = 7 | \M, = dy a A dx a - aA ab dx a Adx b -a { 9 A) ab dy a A dy b , (104) 
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where 9 A is obtained from A by the g-raising of indices, 

( 9 A) ab = g ac g bd A cd . (105) 

We avoid to write simply [A ab ~\ because this might be confused with the 
contravariant inverse of L4 a t], non-existing in our typically physical situations. 
In such situations we have just the following form of (|95|) : 



7||M = 7 ||M S = dy a Adx a . (106) 

And (l98l) becomes then 



n = j (VT y a y b + ~<r ab x a x b + ~ [( 9 a) b a - a) b \ x a y b , (107) 
where the label "<?" refers to the ^-raising of indices, 

(V) Qb := g ac <J cd g db , ( 9 a) b a :=g bc a ca , {»a) a b := a ac g cb . (108) 

The upper-case indices at g refer to the "contravariant inverse", g ac g c b = 5 a t,. 
We write ( 9 a) ab instead of a ah because the latter might be confused with the 
"contravariant inverse" of a ah which is obviously something else than the g- 
shifted object appearing in (|108|) . 

Remark: The first term in (|107p refers to the kinetic energy of the ^-oscillator, 
the second one is its potential energy. The third term is more peculiar. It 
corresponds to something which formally looks like the constant magnetic field 
of induction tensor proportional to 

fab = ja ab = -f ba (109) 

and the covector potential proportional to the linear field 

a k =-fa k ix l . (110) 
In this language the Legendre transformation (|100p has the following form: 

u a = aA ab x b + ^g ab y b , v a = -^g ab x b + aA ab y b , (111) 

and obviously (|101|) (A is eliminated) becomes 

u a = ^gaby b , Va = ~9abX l '■ (112) 

1.4.2 Canonical formalism with "direct" nonlinearity 

Quite independently on our quantum motivation (sometimes more or less hid- 
den), the above study is an interesting and instructive example of how the Dirac 
procedure of Lagrangian constraints works in a rather non- typical situation. One 
might object here against using the heavy formalism for the geometric discussion 
of something that is technically simple: non-dissipative linear finite-dimensional 
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object, formally equivalent to some system of coupled harmonic oscillators, per- 
haps parametrically excited when x does depend explicitly on time. It is well 
known from the theory of linear differential equations in R fc -spaces that techni- 
cally everything is reducible to matrix exponents. Nevertheless, the structures 
revealed and discussed here may be very useful in the infinite-dimensional case. 
And besides, they may be a good starting point towards discussing generalized 
models admitting nonlinearities and second-order time derivatives. Nonlineari- 
ties introduced by the T^-terms in (|4l))) , in particular the ones of the form (|4"9"]) , 
(j50)) , do not modify the above canonical formalism in an essential way. Namely, 
the Legendre transformation is given by the same formulae (|74[) resulting in 
the same primary constraints M (|75|) . ((75)) . According to the formula (fTTj) . the 
expression for the "energy" (|78|) is modified by the V-teivtv: 



e = 7X ab$ a Tp h + V(iP,i?). (113) 
Therefore, the background Hamiltonian ({52")) in (jST)) is replaced by 

n = yxabT^ b + v(ip,iP). (ii4) 



This implies that the vector field (J84J) is modified by the following additive 
correction term: 

dV d dV d 
For example, for the quartic model 

^ = A(r^V) 2 , (H6) 

where A £ R is some constant, we have 

^ = 2A{fT- cd i, d )fT- ba , ^= = 2A (fT- cd i> d ) IW- (H7) 

In virtue of the assumed non-singularity of T, the corresponding equations 
(|8Tj) . (|86|) are again uniquely solvable with respect to the Lagrange multipliers 
at any point of the primary constraints M . The only difference in comparison 
with ([87)1 is that some additional ^-dependent terms appear: 

A« = -pr«W-fr^, r^V^ r - + A|K r - (118) 

2 a 2a g^p° 2 a 2a oip c 

The existence of these unique solutions all over M implies that again the sec- 
ondary constraints M s are identical with the primary ones M, and that they 
are second-class constraints, i.e., M = M s is a symplectic manifold in the sense 
of the structure induced by the symplectic form of the original phase space. 
Referring to formulae ()3"2")) we can simply rewrite (|118|) as follows: 

2 a 2adip a 2 a b 2a g^' y ' 
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with the suggestive, although not very correct, abbreviation: 
dV - dV dV dV - 

|f:=r-^, |L := |Lp* (120) 

8lp a Q^ c dlp^ 9 ^ 

Remark: Let us be careful: as it is mentioned many times, the " Hamiltonian" 
(I114j) is something completely else than the "true" quantum- mechanical Hamil- 
tonian appearing in the Schrodinger equation, even if the non-quadratic term 
V responsible for the direct nonlinearity does not occur at all. Nevertheless, it 
is still responsible for the directly nonlinear quantum evolution of ip as derived 
from our Lagrangian model. 

Using the M-reduced canonical momenta ([90]) 

II a = 2ia$T la , Tla = -2iaT m ip b (121) 

conjugate respectively to ip a and ip and then eliminating the redundancy, we 
find again that according to ([89]) ip and ip are respectively (of course up to nor- 
malisation) the generalized complex coordinates and their conjugate momenta 
on M = M s , namely 

{il> a ,4> h } M =0, =0, {V a ,/} =2^ P *- (122) 

Writing the Hamilton equations of motion on M = M s in the following form: 

dip a 
~dt 

or, equivalently, 

dip 



-ir = {r,H} M (123) 



with the effective Hamiltonian H on M = M s given by ([75]). (1521) and 

making use of all standard properties of Poisson brackets, we obtain respectively 

if^ = H-rf + ir*™, (125) 
at z dip 

-4 " (126 , 

which are evidently complex conjugates of each other. Obviously, just as previ- 
ously, we have put 

a = h, 7 = 2, (127) 

even if our interpretation is purely classical. Those are finite-level Schrodinger 
equations with possibly nonlinear terms controlled by V (ip, ip) . We use the 
term "direct nonlinearity" to stress the fact that T is fixed and the possible 
nonlinearity is just introduced as the perturbation term built of V. 

This is the model suggested by analytical mechanics in the W-space endowed 
with r-geometry of the Hermitian type. Before we go any further towards 
some non-direct, geometry-based nonlinearity of non-perturbative type, some 
comments are necessary. 
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1.4.3 Removing Dirac constraints by second-order terms 

Traditional analytical mechanics of classical oscillatory systems with a linear 
background suggests Lagrangians of the form (j3"5|) , quadratic in tp and perhaps 
corrected by some anharmonic perturbation (the quartic model of the correction 
to the Lagrangian seems to be the most popular one, although, of course, some 
trigonometric, hyperbolic and logarithmic corrections are also used and well 
motivated). On the other side, the Schrodinger-type Lagrangian (|26[) describes 
the usual linear quantum mechanics. Lagrangian (|38p quadratic in derivatives 
does not imply any Dirac singularity; the Legendre transformation is invertible, 
assuming of course (what we do) that F is not degenerate (in standard theory it 
is just positively definite). Therefore, the natural temptation appears to admit 
Lagrangians L(l, 2) (|35p perhaps with a possible anharmonic correction V, when 
it becomes L(l, 2, dnl) (|4"7|) . i.e., the second-order polynomial of velocities with a 
direct nonlincarity. From the Dirac point of view the term quadratic in velocities 
introduces some kind of " regularization" because Legendre transformation is 
then invertible. It is invertible for any, even very small, non-vanishing value of 
|/3|. Without the V-teim one obtains from the variation of ip the second-order 
Schrodinger equation (|4T|) ; the term " Schrodinger" becomes literally true when 
we put (I127j) . i.e., a = H, 7 = 2. The variation with respect to ip itself leads 
to the complex-conjugate equation for ip. In a sense (3 is an additional "Planck 
constant". Admitting a direct nonlinearity via the "potential" V in (|47p. we 
obtain the second-order nonlinear " Schrodinger equation" 

and obviously its complex conjugate 

df Pdrf -6 IdV m 

~ lh ^T ~2^-^ + 23^ r • (129) 
In other words, we have the following expression for the variational derivative: 

*/(l,2,dnl) d^ b dV b 8V 

-= = 2ihT-ab— - pr„ b — - - ^H m i) -=, 130 

5iP {t) dt dt 

where 7 = 2 if for f3 = we are to obtain the usual Schrodinger equation. 
Obviously, performing the variation with respect to ip a , we obtain the complex- 
conjugate expression 

2, dnl) 9 ..# b r R drf -l dV 

sr(t) = ' 2lh ^r r ~ ba ~ p ^ r ~ ba - lijHb ^w (131) 

1.4.4 Some physical interpretation of second derivatives 

It is interesting to mention that there is some motivation from other physical 
problems to admit the term with -^-accelerations, i.e., with the second-order time 



24 



derivatives introduced to the Schrodinger equation. For instance, similar ideas 
were proposed and studied quite a long time ago by A. Barut and more recently 
have been re-investigated by V. V. Dvoeglazov, S. Kruglov, J. P. Vigier and 
others (see, e.g., [TTJIH] and references therein). Among others there is also an 
interesting article where the authors used the analogy between the Schrodinger 
and Fourier equations for nanoscience [5D1 [55] . The short description of the idea 
is presented below. 

Hence, the quantum Fourier equation which describes the heat (mass) dif- 
fusion on the atomic level has the following form: 

= — V 2 T. (132) 



dt 



If we make the substitutions t — > it/2 and T — > tp, then we end up with the free 
Schrodinger equation: 



< i33 » 

The complete Schrodinger equation with the potential term V after the reverse 
substitutions t — > — lit and ifr — > T gives us the parabolic quantum Fokker- 
Planck equation, which describes the quantum heat transport for At > r, where 
t = h/ma 2 c 2 ~ 10~ 17 sec and cr ~ 1 nm, i.e., 

d ^ = ^T- 2 ^T. (134) 
at m n 

For ultrashort time processes when At < r one obtains the generalized quantum 
hyperbolic heat transport equation 



d 2 T dT h n 2V 

— 7T H = — V T ' 

dt 2 dt m h 



r—r -r — = — V 2 T - -£-T (135) 



(its structure and solutions for ultrashort thermal processes were investigated in 
|20j ) which leads us to the following second-order modified Schrodinger equation: 

where the additional term with the second-order time derivative describes the 
interaction of electrons with surrounding space-time filled with virtual positron- 
electron pairs. It is easy to see that (|136[) is analogous to (|4ip if we suppose 
that 

a = h, /3= -4tH, 7 = 2. (137) 

Remark: obviously, an important question appears: is for such a second-order 
Schrodinger equation something like the probabilistic interpretation still possi- 
ble? And if not, what is it to be replaced by? Those are still open questions, 
although some comments will be given below. In a sense the problem is like 
one of the quantum-mechanical interpretation of the Klein-Gordon equation 
versus the Schrodinger or Dirac equations. In our opinion it is impossible to 
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answer all questions at once, so it is a reasonable way to discuss consequently 
any geometrically interesting ideas without waiting for their immediate, per- 
haps premature, physical interpretation. Otherwise one repeats some mistake 
made by Schrodinger who rejected his second-order partial differential equation 
known today as the Klein-Gordon equation. 



1.4.5 Regular Legendre transformation and canonical formalism 

If in (|35|) . (|47p one admits the term quadratic in velocities, i.e., non-vanishing 
0, then the Legendre transformation becomes as follows: 

9L(l,2,dnl) 9L(1,2) .-5 J 

- W> T Ea + 0ip T la , (138) 



di\) a dip a 

-iaTa b ^ b + l3Ta b ^ b . (139) 



9L(l,2,dnl) _ dL(l,2) _ ^ h }b 



dip dip 

The corresponding " energy" function of Lagrangian model is given by 



e = ft + TXaV'V + V (ip, i>) . (140) 



No constraints in the classical phase space appear, so the Legendre transforma- 
tion is invertible and its inverse has the following form: 



r = f«%+ 7 r, ^= r ^-^. (i4i) 

The corresponding " Hamiltonian" is globally defined on the phase space: 



H = ~ (r^aWj; + ia [tt^ - tt^ ) + (jT^ + y X ab^ V- V + V ty, V) . 

(142) 

For any 0^0, this is a regular Hamiltonian and equations of motion may be 
written in the usual canonical form: 

^ Ulfl -h/X d7Ta I*r H\ C\A.*\ 

~dT = {1P ' H} = ~dt = - ~W ( 5 

or in the equivalent form of complex-conjugate quantities: 



I J dn-a at Q^ a 



There are no phase-space constraints and no Lagrange multipliers. 

The limit transition — > is a rather complicated and obscure problem. The 
equations of motion (|128[) , (|129[) transform then smoothly into (|125[) . Q126p . 
Obviously, this concerns the very form of equations, but the phase portraits 
change catastrophically as it is usually the case when the highest-order derivative 
terms in equations are neglected. The energy expression (|140[) also reduces 
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smoothly to (|113[) . But it is no longer the case with the Hamiltonian (|142[) . 
There is no directly well-defined limit when (3 in (| 142|) tends to zero. In any 
case nothing like (f5Tj) , {82]), (jTT5J) appears. This is due to the fact that for 
f3 = the Legendre transformation is not invertible. For f3 ^ the number of 
complex degrees of freedom remains n, i.e., that of real ones remains 2n. But 
for the vanishing j3, the number of effective real degrees of freedom becomes n. 
This qualitative discontinuity resembles similar phenomena in field theory, e.g., 
passing over from the non- vanishing mass of " photons" in the Proca equations 
to the vanishing one in the Maxwell ones. 

In this sense the non- vanishing f3 leads to some " regularization" of the un- 
derlying Hamilton mechanics, but it is not yet clear how this regularization 
interferes with the probabilistic interpretation, the structure of the effective 
scalar product of "wave functions" and its positive definiteness. This resembles 
the corresponding problem for the Klein-Gordon equation. 

2 Non-direct nonlinearity and Hamiltonian sys- 
tems on manifolds of second-order tensors 

The mentioned above models of nonlinearity in a finite-level Schrodinger equa- 
tion or second-order " Schrodinger equation" were " direct" in the sense that the 
term of Lagrangian responsible for nonlinear phenomena was a kind of perturba- 
tion extra introduced as an additive correction to the linear background. There 
is a rather large freedom of a priori admissible models and no very convincing 
criteria of choice do exist. The choice is indeed a matter of intuition and is based 
on a kind of phenomenology. And because of this it is never very convincing. 

Below we try to discuss some models where nonlinearity is not introduced 
"by hand", but rather it is based on some kind of geometric aprioric arguments. 

2.1 Removing absolute objects and admitting dynamical 
scalar product 

Qualitatively, the idea is as follows: Lagrangians (|26|) . (j35l) . (f36|) . (|46|) . (|47)) 
contain only one dynamical quantity tp with n complex degrees of freedom, i.e., 
2n real ones (x,y). As we saw, it is convenient to follow the method popularly 

used in field theory and consider formally the components ip a and their tp as 
independent variables. But all the mentioned Lagrangians contain also some 
absolute (non-dynamical) object, namely, the scalar product V . It is fixed once 
for all, just like the metric tensor of Euclidean space or that of Minkowskian 
space-time of special relativity. And it is just here some doubts appear. Namely, 

Nature does not like absolute objects. 

It turned out that in general relativity the metric tensor of our four-dimensional 
(in general curved) space-time is a dynamical quantity which together with 
"physical" or "matter" fields satisfies a closed system of differential equations. 
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This system is based on the mutual interaction and the metric tensor turns out 
to be also some kind of a physical field, namely, it just describes relativistic 
gravitation. Hence, there are no absolute objects in space-time. 

One can ask why not to follow this pattern and admit T to be a dynamical 
object mutually interacting with ip. For what would be the physical reality hidden 
behind some absolutely fixed T ? 

It turns out that the only natural Lagrangian dynamics of the T-object is 
strongly nonlinear and this introduces a kind of effective nonlinearity to the 
dynamics of the total (T, ■0)-system. This nonlinearity has geometric origin, 
therefore, it is well motivated and non-perturbative. By that we mean that it 
is not a correction to any well-defined linear background and because of this 
it cannot be analysed by standard techniques of the perturbation calculus, any 
expansion with respect to some "small" parameter, etc. 

Analytically T is represented by the quadratic matrix with coefficients r s t,. 
Geometrically it is a Hermitian element of the tensor space W (8 W* , i.e., a 
Hermitian form on W . In various problems of analytical mechanics and math- 
ematical physics one often deals with dynamical systems the phase spaces of 
which are byproducts of some matrix manifolds. Obviously, those matrices ap- 
pear on the level of calculus and analytical representation but as a matter of 
fact they represent various second-order tensors, i.e., linear mappings operating 
between two linear spaces. In applications which present interest for us they 
are usually isomorphisms between linear spaces of the same dimension, i.e., we 
deal with quadratic non-singular matrices. Some exceptional models when one 
deals with rectangular (not necessary quadratic) matrices are presented in [25] . 

2.1.1 Manifolds of linear mappings 

Let us begin with some introductory remarks. Matrices provide analytical de- 
scription of linear mappings. So, let U, V be some linear spaces of the same 
dimension. At a moment, we do not precise if they are meant over the real of 
complex field. Let L(J7, V) denote the linear space of all linear mappings of U 
into V and let LI(f7, V) C L(C/, V) denote its open submanifold consisting of 
linear isomorphisms of U onto V. When dealing with analytical formulae, we 
shall use some bases (. . . , Ea, ■ ■ ■), (■ ■ ■ , e*, . . .) respectively in U and V. Linear 
mappings ip G L({7, V) are analytically represented by matrices with coefficients 
ip 1 A meant in the following convention: 

ipE A = e^A- (145) 

The conjugate mappings ip* : V* — > U* is given by 

<p*p = po<p, <p*e i =y i A E A , (146) 

where (. . . , E A , . . .) , (. . . , e\ . . .) denote as usually the dual bases of dual spaces 
U*, V*: 

E A (E B ) = (E A ,E B )=6 A B} e i (e j ) = (e i ,e j ) = S i j . (147) 
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And as usual we make use of abbreviations: 

L(Z7) := L(U, U), L(V) := L(V, V), etc. (148) 
GL(f7) := LI(U, U), GL(V) := LI(V, V), etc. (149) 

Obviously, L({7), L(V) are associative algebras and at the same time commuta- 
tor Lie algebras. They are of course canonically identical with Lie algebras of 
the corresponding Lie groups GL + (U), GL + (V), i.e., the connected components 
of GL(U), GL(V). 

In applications we often deal with various closed submanifolds of L(U,V), 
usually distinguished by some additional structures in U, V like, e.g., scalar 
products. First of all, let us stress that the groups GL({7), GL(V) act in a 
natural way on LI({/, V), respectively on the right and on the left. So, for any 
A £ GL(V), B £ GL{U) we have the action: 

U(U, V)3y^ AipB = {L A o R B ) {if) = (R B o L A ) (tp), (150) 

where the symbols La, Rb refer respectively to the left- and right-hand-side 
actions. Obviously, both Lgl(V)j Rgl(U) ac t effectively and freely on Ll(U, V). 
However, it is important that GL(V^) x GL(i7) does not act effectively through 
^GL(y)^GL((7), because dilatations in GL(V) and GL(U) act in the same way 
on LI([/, V) and the kernel of non-effectiveness of GL(V^) x GL(U) in the action 
through (I150P is given by the subgroup 

{(AM v ,A -1 Idtr) :AeR\{0}}. (151) 

Let us mention that L(f7), L(V) are also semigroups under the composition 
of mappings and these semigroups act on the total L(i7, V) according to the 
formula (1150[) . Obviously, if A, B are not isomorphisms, i.e., elements of GL(V), 
GL(J7), the submanifold Ll(U, V) is not preserved by them. If linear spaces 
U, V are endowed with some scalar products, i.e., symmetric (usually non- 
degenerate) bilinear forms, -q € U* ®U* , g £ V* ®V* , then the special attention 
is paid to submanifolds of "isometries" 0(£7, rj; V,g) C Ll(U,V), consisting of 
such mappings ip £ Ll(U,V) that 

9 = <P*V, 9tj = VAB(fi A iip B j- (152) 

They are homogeneous spaces of isometry groups 0(V, g), 0(U,rj) acting in the 
sense of (1150[) . These isometry groups consist of transformations preserving 
respectively g and 77, 

g = A*g, V = B* V , (153) 

i.e., analytically, 

9ij = g k iA k t A 1 ^ r, KL = vmnB m k B n l . (154) 

In such applications we usually deal with real linear spaces and symmetric pos- 
itively definite metrics. Nevertheless, everything works well at this stage also 



29 



for general tensors 77, g, even without any kind of symmetry, or, in the extreme 
case, when 77, g are skew-symmetric, i.e., if one deals with symplectic structures 
in U, V. Another interesting case is when one considers as relevant subman- 
ifolds of U(U,V) the orbits of the subgroups U(V) x U(£f), SL(V) x SL(U) 
acting through ((T5D|> . Here U(V), U(U) and SL(V), SL(U) denote respectively 
the unimodular and special linear subgroups of GL(V), GL(U), i.e., groups con- 
sisting of linear mappings the determinants of which respectively have modulus 
one or are just themselves equal to one. If U, V are real linear spaces, then 
U(V), U(f7) preserve the volumes, i.e., the Lebesgue measures in V, U (it does 
not matter how normalised), and SL(F), SL(U) preserve both the volumes and 
orientations (both standards of orientation separately). The manifold LI(U,V) 
is then foliated into orbits which are homogeneous spaces of the mentioned sub- 
group acting through (|150|) . In the case of real spaces the strata of the action of 
unimodular groups have two connected components corresponding to the fate 
of orientation under the group action. 

2.1.2 Hermitian metrics and unitary mappings 

If linear spaces U , V are complex, then in applications we are interested rather 
in sesquilinear Hermitian products 77, g than in bilinear ones. The orthogonal 
groups are then replaced by the unitary ones U(f7, 77), U(V,g). They consist of 
transformations A £ V(V,g) C GL(V), B G U(U,rf) C GL(U) satisfying ((155)1 
in the analytical form: 

% = 9kA k i Al h nwL = nMN~B M KB N L. (155) 

The Hermitian metrics g G V * <S> V* , 77 G U <£> U* give rise to the manifold of 
unitary isometries U(C7, 77; V, g); it consists of linear mappings relating 77 to g, 

V = <P*9, njB = ffij^A^B- (156) 

Obviously, the subset U(Z7, 77; V, g) is a homogeneous space of both U(V, g), 
U(U, 77) acting through (|150j) . In applications we are interested in, one deals 
usually with Hermitian positively definite metrics 77, g, but some general state- 
ments are valid without this restriction, i.e., pseudo-unitary groups are formally 
admissible. Similarly, if U, V are real linear spaces, then usually 77, g are sym- 
metric positively definite forms. It is so in mechanics of affinely-rigid bodies 
investigated by us and others in various papers (5) (6) (7) [TBI EH 1221 O (3U [35] 
(then U , V are translation spaces respectively in the material and physical affine 
spaces) . 

2.2 Non-holonomic velocities and invariants 

Let Q C LI(f/, V) be a submanifold used as a configuration space of some 
"analytical mechanics". If U, V are complex, then Q may (however need not 



30 



to) be analytic. Moreover, in applications we have in mind, it usually is non- 
analytic. "Motions" are described by curves 

R9 *-►¥>(*) eQcU(U,V). (157) 

Both from the point of view of geometrical foundations and practical calcula- 
tions, it is often convenient to use non-holonomic velocities which at the time 
instant t € K are given by 

n(t) := ^(^(ijGLfy), (158) 
n(t) := ^- 1 (*)^(*) = ^- 1 (*)n(*M*)eL(tT), (159) 

(XI 



i.e., analytically, 



n 3 = —^-(p j, n B = ip i—jj-- ( 16 °) 



If the manifold Q C LI(C/, V) is parameterized by generalized coordinates q^, 
11 = 1, ... , dim Q, then obviously 



0' i ( 9) «) = fiV(g)^, ^bM) = ^(<z)^ (161) 



where 



i d(p l A ^-XA 5A ,-lA 9(f'B nR0 , 



The advantage of using the quantities il, fl instead of the generalized velocities 
is that due to their tensorial structure one can construct of them some scalar 
invariants: 

I p := Tr (n p ) = Tr , p = 1, . . . , dim V; (163) 

according to the Cayley-Hamilton theorem, taking other values of p one does 
not obtain anything new, just some functions of (|163p . Obviously, the most 
important are expressions quadratic in velocities, i.e., Ii and , because in 
analytical mechanics they are used for constructing kinetic energy models {Ii is 
the main term and I\ is a merely correction). Even for some non-quadratic (in 
velocities) hypothetical models of the kinetic energy T, the quantities I p , first 
of all I2 and (ii) 2 , are reasonable modules for constructing T in a form, e.g., 

T = f (l 2 ,{hf) , (164) 

where / is some appropriately postulated function of two variables (incidentally, 
let us mention that in the three-dimensional spatial formulation of relativistic 
mechanics of the material point the kinetic Lagrangian is an irrational function 
of the squared absolute value of velocity). 
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2.2.1 Kinetic energies, i.e., Riemannian metrics on manifolds of lin- 
ear mappings. High-symmetry models 

Quadratic forms of velocities tangent to Q C LI(£7, V) are geometrically equiv- 
alent to some metric tensors on Q. Usually, at least in our problems, Q is 
somehow special from the point of view of geometry of LI([7, V) and the most 
natural metrics on Q (viable models of kinetic energy) are restrictions to Q of 
geometrically distinguished metric tensors on the total Ll(U,V). And this is 
where we must begin with some digression concerning our former work on an- 
alytical mechanics of afhnely-rigid bodies, because some ideas developed there 
seem to be an inspiration for our search of geometric Schrodinger nonlinear 
models. Let us remind that if U, V are real linear spaces of material and spa- 
tial translations respectively, then the usual formula for the kinetic energy of 
internal (relative) degrees of freedom reads [3S[ [33J 

1 d^A^BjAB (165) 

2 yu dt dt y ' 

where g £ V* ® V* is the metric tensor of the physical space and J £ U ® 
U is, roughly speaking, the constant co-moving tensor of inertia (the second- 
order tensor moment of the mass distribution with respect to the Lagrange 
coordinates). This expression is invariant under (|150p . where A, B are confined 

respectively to 0(V,g), Ou7, jj. In the special case of inertially isotropic 

affine body, when J is proportional to the reference material metric rj in the 
contravariant reciprocal form, thus, 

J AB =Irf B , J AC J CB = 6 A B , (166) 

the expression (|165p is isotropic both spatially and materially. If ip is an isom- 
etry, i.e., ip £ O (U, rj; V, g), then (|165[) becomes the usual kinetic energy of 
the metrically-rigid body (gyroscope) and the resulting metric tensor is just 
the restriction of the metric underlying (|165[) in LI(£7, V) to the submanifold of 
isometries. Just as the latter one, the metric of (|165l) is not invariant under the 
total group (|150p . And it is just here where the quasivelocities (|158l) - (|160[) be- 
come interesting. Namely, they transform under (|150|) according to the following 
rule: 

Q^ACIA- 1 , h-tB^hB. (167) 

Because of this, the corresponding quantities I p (|163p are invariant under the 
total (|150p . i.e., under the most natural group underlying geometry of degrees 
of freedom. Therefore, the corresponding kinetic Lagrangians built of I p , in 
particular the quadratic ones built of I2 and (Ii) 2 , and their underlying metric 
tensors are also affinely invariant, i.e., non-sensible to the action of (|150p . The 
underlying metric tensors on LI([/, V) are essentially Riemannian (they have 
non- vanishing Ricmann tensors). This curvature is due to the fact that in ex- 
pression for T given by combinations of I2 and with constant coefficients, 
the corresponding quadratic forms of generalized velocities dip 1 a 1 1 dt have irre- 
ducibly (^-dependent coefficients. And no change of generalized coordinates for 
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some new ones may help here. The corresponding quadratic forms of if 1 will 
have always g-dependent coefficients. The reason is that ft, ft are essentially 
non-holonomic quantities and this in turn follows from the non-commutativity 
of linear groups GL(V), GL(U). 

Proposition: As mentioned, from the point of view of geometric a priori, the 
most natural models of kinetic energies have the following form: 

T = ±Tr (ft 2 ) + f (Tr ft) 2 = ^Tr (ft 2 ) + f (lr ft)* . (168) 

This expression and the underlying metric tensor on IA(U, V) are invariant under 
the total (|150[) . And this is the most general metric of this property. This is an 
affine counterpart for the kinetic energy of the spherical metrically-rigid body. 
Proposition: Of course, in principle, just like in mentioned rigid body me- 
chanics, one can think about models affinely invariant in space or ones affincly 
invariant in the material of the body. As in 32, 33] they are respectively given 

by 

T = ^£ B A D c n A B h c D , (169) 
T = l-RPh&j&u (170) 
obviously with constant coefficients C, 1Z. 

Remark: The only situation of (|169[) and (|170[) to coincide, i.e., of the in- 
variance under the total (|150[) . is just (|168p , the geometrically most natural 
situation. One should stress that (|168[) is not positively definite, its main term 
(A/2)Tr (ft 2 ) has the signature (n(n + l)/2 +,n(n— l)/2 — ), so, if gyroscopic 
constraints are taken into account, A must be negative. It was shown in our 
earlier papers on affine bodies that the non-definiteness of (|168[) may be just 
convenient in certain models of elastic vibrations; in those models the dynam- 
ics of deformative motion is encoded in the kinetic energy term, i.e., in some 
kind of the effective metric on LI({7, V). This resembles in a sense the idea of 
Maupertuis principle. 

Obviously, in certain phenomenological models of deformative dynamics and 
the rotation-deformation coupling it may be reasonable to try to mix (|168p (or 
(|169p and (|170[1 ) with the metrical term Q165p , especially with its isotropic form 

Proposition: In certain problems it may be perhaps reasonable to postulate 
the affine model f|168[) perturbed by two corrections breaking the background 
affine symmetry to the orthogonal one, both in (V,g) and (U,rf), e.g., [3"2l [55] 

T = ^jW'nVfi*, + jmLV MN ^ K M^ L N + ^Tr (ft 2 ) + I (Tr ft) 2 (171) 
= yW'flW + jVKLV MN ^ K Mn L N + ^Tr (ft 2 ) + | (Tr ft) ' . 
If we use orthonormal coordinates in which gtk =* Sik, t\ab =* Sab, then 
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the following simple formula based on the R ra -object may be used for (jTTTJ) : 



T = ^-Tr (n T n) + j-Tr (fi T fi) + |Tr (Q 2 ) + | (Tr tt) 2 

= ^Tr(r! T r!) + |Tr(^) +^Tr(^ 2 ) +|(Tr^) 2 . (172) 

Remark: Let us stress the following important circumstance. Unlike (I17ip . 
(1172j) . the expression ()168|) does not preassume any fixed metrics in U, V. This 
will be just the pattern to be followed in our model of "non-direct" geomet- 
ric nonlinearity in dynamical systems motivated by the Schrodinger equation. 
There we deal, of course, with complex linear spaces. The above formulae, 
although providing some guiding hints, cannot be literally used because the 
corresponding Lagrangians would be either trivial or complex. Neither their 
real nor imaginary parts 

ReT=i(T + T), Imr=^(T-T) (173) 

would be useful; they do not correspond to any expressions interpretable in 
quantum-mechanical terms. 

As yet, the linear spaces U, V were assumed completely unrelated to each 
other. And it was just correct in the mentioned applications to clastodynamics. 
Let us remind that L(£/), L(l/) are canonically isomorphic with the commutator- 
sense Lie algebras of GL(U), GL(V) and the expressions (|163j) were their Casimir 
invariants; in particular, the appropriate special case of (I168p . i.e., 

A = 2n, B = -2, (174) 

corresponds to the Killing metric (degenerate on the total L([/), L(V) because 
those algebras are not semisimple; dilatations form the normal subgroups of 
GL(t/), GL(V)). In the mentioned elastodynamical applications we often deal 
with the situation where the configuration space Q C LI (U, V) is an orbit of 
some subgroups G v C GL(C7), G v C GL{V) acting through (fT5U|l . Those 
subgroups are isomorphic and for any <f) S Q we have that 

Gu - ^Gycb, G v = cbGucjj- 1 . (175) 

Remark: The usual, i.e., metrically- rigid, body is a typical example, we have 
then 

Gu = 0(11,7]), G v = 0(V,g), 4>eO(U, m V,g). (176) 

More precisely, in some realistic mechanical applications we are dealing then 
with the connected components 0+ (U, -q) = SO(C/,77), + (V,g) = SO(V,g) 
and their orbits. Similarly, one considers orbits of SL(£7), SL(U) (incompressible 
body), etc. Restricting the expressions (|163p to Lie subalgebras G' (J C L(Z7), 
G' v C L(U) one obtains some Casimir invariants; usually they need not be 
independent. For example, for SO (U, 77), SO (V, g) we have I p = for any odd 
P- 
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Before going any further we must discuss some special situations. As men- 
tioned, the linear spaces U, V above where independent on each other. The 
corresponding configuration spaces consisted of linear mappings ip £ ~L(U,V), 
i.e., of tensor quantities <p G V ® U*\ they were analytically represented by ma- 
trices <~p % a (doubled quantities in the Schouten-Veblen language). But when dis- 
cussing Hamiltonian systems inspired by the Schrodinger equation, we must use 
matrices which analytically represent some scalar products, i.e., twice covariant 
tensors in W, more precisely, the sesquilinear Hermitian forms T € W <g> W* 
represented analytically by matrices [Tab] (in realistic applications positively 
definite ones). Some relatively new structures appear then. 

2.2.2 Metrics on groups as the special case 

Before discussing the manifolds of scalar products and the aforementioned struc- 
tures, it is instructive to start, for comparison and for preparing the proper 
mathematical instruments, with the situation U — V. Then L1(U, V) becomes 
simply the linear group GL(V) — GL(U) and the resulting scheme is that of 
(more or less) invariant Hamiltonian systems on the total Lie groups or their 
subgroups with some interplay of left and right invariance. Configuration space 
consists of non-degenerate mixed second-order tensors (once contravariant and 
once covariant) in a given linear space. The corresponding metric tensors under- 
lying the kinetic energy expressions like (|169|) , (|170[) and their special cases like 
(|168|) . (I171|) have respectively the following explicitly non-Euclidean (if n > 1) 
forms: 



Q = £ j n l m v- ln i v- lm k d l p i j ®d V > k l , (177) 

g = y- lj n y- ll m n n i m k d V i j ®dy k u (178) 

G = {Ap- u l p^ k +Bp- l3 l p- ll k )dp l ] ®dp k u (179) 

g = htp^m^ng^gikdip'j <g> d(p k i + l2gmn^ lm i ^' ln k g jl dip i J (g> dip k i 

+ (Aip- U i<p-V k + Bip- l U<f~ ll k) dtp'j ® dip k h (180) 



Similar expansions with respect to the basic terms dip 1 a <8> dip 3 b m ay be done 
in general, e.g., for (|171[) when U, V may be different linear spaces. The corre- 
sponding formulae are structurally like the above ones. The explicit expressions 
like (|177p - (|180p are in spite of their rather technical nature interesting in them- 
selves and give an alternative analytical insight into the structure of expressions. 

The restriction of the above metrics/kinetic energies to submanifolds of 
GL(V) is analytically achieved by specifying ip 1 j as functions of some parame- 
ters, i.e., generalized coordinates q 1 *, /i = 1, . . . , /. Technically this implies that 
the differentials dpi 1 j are specified as 

- (181) 

In applications the mentioned submanifolds of Gh{V) are usually its subgroups, 
e.g., special orthogonal SO(V, g), special linear SL(T^), etc. In formulae one uses 
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Lie algebras of the mentioned subgroups as geometrically and technically con- 
venient non-holonomic velocities. Usually the choice of some particular mod- 
els of T, i.e., of Q, is motivated by the particular symmetry demands under 
left and right regular translations in GL(V) or in the corresponding subgroup 
G C GL(y). The highest symmetry corresponds to (|179|) . 

2.2.3 Manifolds of scalar products and their non-holonomic veloci- 
ties 

Let us now consider another pair of mutually related spaces, namely, the dual 
pair; V = U*, U = V* . Linear mappings from V to V* are twice covariant 
tensors a, i.e., elements of V* ® V* ~ L (V, V*), represented by matrices ay. 
Similarly, linear mappings from V* to V are twice contravariant tensors (3, i.e., 
elements of V ® V — L(V*,V), analytically represented by matrices (3 13 . In 
the case of isomorphisms there exists a natural canonical bijection of LI (V, V*) 
onto LI (V* , V); namely, according to the standard rules, interrelated are objects 
a e V* <g> V*, (3 S V <g> V such that 

a lk ^=8 % \ f3 ik a kj =S i j . (182) 

One can write simply 

P = oT 1 , a = (3~ 1 . (183) 

Just as when dealing with geodetic Hamiltonian systems on GL(V), one consid- 
ers now again the very special case of LI(f, V), however, with completely new 
physical and geometrical peculiarities. 
For any time evolutions 

l3in a(t) e LI (V, V*) , Rsfn (3{t) e LI (V*, V) (184) 

the fl, fi-objects are as well defined as for any general R 3 t h- > tp(t) € LI (U, V). 
However, certain new features and peculiarities appear. The corresponding 
objects for curves in LI (V, V*) are denoted by fl[a] € L (V*) ~ V* ® V and 
f2[a] € L (V) ~ V®V*; more explicitly one can write them as (a, a), (a, a). 
Analytically they are given by 

n[a\i j = a ik a- lk i, h[a] 1 = cT uk 6t k] . (185) 

Obviously, they are interrelated by the following formulae: 

n[a\i j = a ik n[a] k ia- llj , Q[a] = an[a]a -1 , (186) 
n[a]%- = a- lik n[a] k l aij, fi[a] = a -1 n[a]a. (187) 

Similarly, for curves in LI (V*,V) we use the symbols Cl\J3] € L (V) ~ V" ® V*, 
f2[/9] € L (V*) ~ V* ® V or (p, ^ , . The corresponding analytical 

expressions are as follows: 

nffij = $ lk p-\ 3 , - p-\J kj . (188) 
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In analogy to (|186p , (|187[) there exists the following obvious relationships: 



(189) 
(190) 



One can easily show that 



n [/r 1 ] 



n a 



1 



(191) 
(192) 



We are dealing here with the special case of the general scheme of LI (U, V)- 
models, therefore, the properties of a £ V* <g> V* ~ L (V", V*) or /3 g V ® V ~ 
L (V*, V) may be analysed from the point of view of some additional structures 
like metrics in V, V*, in analogy to metrics j £ V* ®V* , i] G £/* <g) [/* in 
mechanics of affine bodies. In particular, one can consider the rigid motion in 
the sense of those metrics; the Q, £!-objects become then skew-symmetric with 
respect to introduced metrics, i.e., are interpretable as "angular velocities". 
Obviously, when dealing with the manifold of scalar products in V, it would be 
rather artificial and exotic to introduce two independent metrics in V and V* 
(as we did in mechanics of affine bodies in 1A(U,V)); rather some metric in V 
and its contravariant inverse in V* would be used. 

Remark: One point is important here. In mechanical theory of systems with 
affine degrees of freedom, in LI(U,V), GL(V) as configuration spaces, it was 
rather natural to discuss constrained motion along subgroups of GL (V) or some 
orbits of the left or right actions of subgroups of GL (V), GL (U) on LI (U, V). 
Nonholonomic velocities were then the elements of the corresponding Lie sub- 
algebras. And one concentrated on the left or right (or both) invariant metrics 
(kinetic energies) on configuration submanifolds. Nothing like this is useful in 
applications of dynamical systems on the manifold of scalar products. For bi- 
linear forms, i.e., elements of V* <g> V* or V (8 V, the symmetry /antisymmetry 
is well defined without any reference to something like a once fixed absolute 
scalar product. Usually we deal with symmetric or antisymmetric scalar prod- 
ucts, e.g., if V is over reals, with (pseudo-) Euclidean spaces or (generally) with 
symplectic spaces. When V is over the complex field C, the special stress is 
laid on Hermitian scalar products. Obviously, being twice covariant (or twice 
contravariant) tensors, not mixed ones, such objects cannot be multiplied in 
spite of their analytical matrix form. Even if we identify V with V* using some 
pre-fixed reference scalar product and so identify the mentioned forms with lin- 
ear transformations, nothing like the subgroup structure survives because, as a 
rule, the subsets of symmetric, antisymmetric, or Hermitian matrices are not 
closed under multiplication. And as a rule, the above objects O, f2 do not form 
Lie algebras. Nevertheless, they are well-defined mixed tensors in V and enable 
one to construct invariant quadratic scalars and, more generally, the family of 
basic scalars homogeneous in derivatives (generalised velocities) in a complete 
analogy to Lie-algebraic Casimir invariants (|163|1 . 
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2.2.4 Canonical Riemann structures, i.e., kinetic energies on mani- 
folds of scalar products. High-symmetry models 

So, we would like to fix now some useful (at least hopefully) models of "ki- 
netic energy", i.e., Riemannian structure, on appropriate manifolds of scalar 
products. As usual, it is impossible to be at the same time very general and 
computationally effective. And as a rule, it is the special cases, first of all the 
ones with high symmetry, that has a chance to be physically viable. So, step by 
step one reduces the interest to (|169[) . (|170l) . later on to (|171[) . (|172p . then to 
(I168j) and first of all to its special case B = (it is clear that the square-term 
controlled by B is a merely secondary correction). 

Explicitly, the Riemannian structure on the manifold Sym (V* <8> V*) C V*® 
V* of symmetric scalar products on V, constructed by analogy with the GL(V)- 
prescription (|168[) . has the form corresponding to the kinetic energy 

T=^Tr(f7H 2 ) + |(Trfi[a]) 2 , (193) 

i.e., analytically to 

T=^{a} i m{a}/ + ^(n{a} l f. (194) 

In analogy to (|179[) and using the standard Riemann expressions, we write for 
the underlying metric tensor that 

Q = {Aa- m a~ 1]k + Ba- ljl a~ llk ) da l3 ® da M , (195) 

notice however some essential differences between (|195l) and (|179p . Using the 
favourite physicists way of thinking, we have the metric element 

ds 2 = (Aa- m a- ljk + Ba- l3i a- llk ) da i3 da ku (196) 

where the pairs (ij), (kl) are, roughly speaking, bi-indices; their ordering does 
not matter when we deal with the manifold Sym (V* ® V*) of symmetric forms 
on V. 

It is clear that when the forms a are symmetric, the plenty of "aesthetic" 
changes of ordering of indices is possible. Something similar, although a bit 
different, may be done for manifolds of symplectic forms. Equation (|195[) may 
be written in the following form: 

g = g» kl (a)da ij ®da kU (197) 

where, let us notice carefully, 

gabed = _ ^ a -lac a -lbd + ^-lbc^-lad) + Ba -lab a ~lcd _ ^gg) 

The above formula implies that g abcd have all necessary symmetry properties to 
represent some Riemannian metric on the manifold of real (or complex-analytic) 
symmetric scalar products, 



iabcd f-rbacd r->abdc r-icdab 



(199) 
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(strictly speaking, the symmetry under the simultaneous exchange of pair (ab), 
(cd) holds without the representation (|198[> ). 

Although not interesting for our purposes here, it is nevertheless interesting 
in itself to search for natural metrics on the manifold of all non-degenerate 
bilinear forms on V, not necessarily symmetric ones. They work as follows 
when evaluated on pairs of tangent vectors: 

G abcd u ab v cd , (200) 

where 

gabcd = Aa -lab a -lcd + tL {oT lac oT lhd + a" 1 "^" 1 ^ 
i ~D « — lad _ — lcb i / —lad _ —lbc , _ — Ida _ —lcb\ 

+ dol a H [a a + a a 



D 






-l&d 




~2 


(a- 


lac cT 


+ a 


E 






-16c 




~2 


(0- 


lo V 


+ a 


F 


(a~ 


lca aT 


■lbd 


+ a' 


~2 









+ Ga' lda a- lbc + -(a- lca a- lbd + a- lac a- ldb ). (201) 

Obviously, in an exactly the same way as we did above, we can construct natural 
scalar products on manifolds of non-degenerate twice contravariant tensors, i.e., 
Riemann structures on Sym (V ® V). Thus, we use the /3-tensors, and, e.g., 
instead of formulae 1(155) ) -1(155 )1 we obtain respectively that 

T=|Tr(f7[/5] 2 ) + |(Trfi[/3]) 2 , (202) 

i.e., analytically 

T = ^n[0] i j n[f3y i + ^(n[i3} i i )\ (203) 

The underlying metric has the following form: 

Q = + Bp£p£) dp* <8 d(i k \ (204) 

in a full analogy to 1(155)1 . And just like in l(T57|) . (fT^g)) we have that 

Q = Q l0 ki{P)d^ ®d!3 k \ (205) 
Gabcd = yfe^+^^+B^ft 1 - (206) 

All these metrics on the manifolds of scalar products are evidently curved, just 
like the Killing tensors on semisimple Lie groups. 

Remark: Let us notice that if some fixed metric G € Sym (V* <8> V*) is distin- 
guished, it does not matter why, then the analogues of (|17ip . (|172p . (| 1 80|) are 
also well defined both on Sym (V* ® V*) and Sym (V ® V). Analytically, the 
counterpart of (|180|) as a metric on Sym (V* ® V*) is given by 

= ha- ljm a- lln G mn G ik d ai j ® da fc ; + l 2 oT Um a~ lkn G mn G fl da^ <g> da kl 
+ (ia" 11 ^- 1 ^ + Ba^J'a-"*) da^ ® (207) 
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where I\, I2, A, B are some constants. 

Obviously, G with the upper-case indices is the contravariant inverse of G 
with the lower-case indices, G lk Gkj — S l j. One must be very careful when the 
inverse symbol is omitted for brevity in a~ l1 ^ , i.e., when we use the simplified 
notation a % K This is again the inverse, 

a ik a kj =6 i j , (208) 

however, one must remember that it is something else than the G-raising and 
lowering of indices, 

a ij = a -Uj ^ G tk G ]l a M . (209) 

If nothing like some distinguished G is fixed, we can consider only geodetic 
models, first of all homogeneous ones, characterized by high symmetries, like 
(|193|) - (|195p . If for some physical reasons some certain reference metric G is 
fixed, we have more possibilities. First of all, we can take the model (|207p of 
the kinetic energy and manipulate somehow with the constants I\, I2, A, B. 
But there are also some natural classes of potentials V(a) and the corresponding 
Lagrangians L = T — V(a). Namely, having at disposal two "metrics" G,a £ 
Sym (V* ® V*) we can construct the mixed tensor a <E V® V* ~ L(V). Then in 
the n-dimensional space V we can construct the system of n independent basic 
scalars: 

I p (a, G) := Tr (a p ) , p=l,...,n. (210) 

According to the Cayley-Hamilton theorem, I p for any other integer p may be 
expressed as a function of the above ones; again the property which was used 
in the study of deformation invariants. And the most natural and symmetric 
potentials V are appropriately chosen functions of I p (a,G). 

Obviously, exactly the same may be done for systems on Sym (V (g) V); there 
is no need to write down the obvious formulae. 

2.2.5 Taking "translations" into account. Riemann structures on 
the manifolds of wave functions times scalar products 

Let us also mention another class of canonical Riemannian structures on the 
manifolds V x Sym (V* ® V*). They have the following form: 

Q = \ !<>, ,< In' O du 3 + (Aa- 1H a- ljk + Ba~ lji a - llk ) da %3 ® da ku (211) 

where M 1 A, B are constants and u l are linear coordinates on V corresponding 
to our choice of basis (therefore, they are simply elements of the dual basis in 
V*). Kinetic energy based on (|21ip is as follows: 

T = ™ o^V + ^Tr (ttH 2 ) + I (Tr Q[a}) 2 , (212) 

where v l = du l jdt is the "translational" velocity, in analogy to formulae for the 
affinely-rigid body. 
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Remark: Let us observe that in the " translational" part of T the velocity 
v l is squared with the use of a itself, not with the use of some fixed metric 
G € Sym (V* ® V*). This resembles some models of affine motion, where the 
translational velocity is squared with the use of the Cauchy deformation tensor 



2.3 Natural Lagrangians of high symmetry, dynamical 
equations for scalar products and simple solutions 

Those were preliminary remarks based on intuitions developed during our earlier 
study of Hamiltonian systems on groups and homogeneous spaces, first of all on 
manifolds of affine and linear mappings (affinely-rigid bodies). Now we are well 
prepared to return to our proper subject, i.e., to the non-direct nonlinearity of 
hypothetical quantum mechanics. 

Let us go back to the complex linear space W and consider the manifold of 
Hermitian scalar products there, i.e., the manifold of non-degenerate sesquilin- 

ear Hermitian forms, Herm (w ® W*J . In principle they should be positively 
definite, but in many problems this restriction is not formally necessary. Ana- 
lytically such scalar products T are represented by Hermitian matrices [Tab] in 
the sense that 

T(u, v) = T (u a e a , v b e b ) = r s6 uV\ (213) 

where e a £ W are some basic vectors in W. 

Let Herm (w* ® W*^j, or rather some its connected component (first of all 
the one consisting of positive forms), be our configuration space. Obviously, 
in spite of the complex character of W, the set Herm (w* (8 W*^j is a real 
linear space, and the mentioned configuration space is a real manifold, an open 
subset of Herm (w*®W*y We are interested in Ricmannian structures on 

this manifold and mainly in ones analogous to (|193[) - (|195[) and (|202[) - (|204p . 
To be more precise, those Riemannian structures will be restrictions of some 
Hermitian ones defined on the total W ® W* . We introduce them as some 
kinetic energy forms. 

Proposition: Obviously, the only natural counterpart of (|193fl — <|195[) is as fol- 
lows: 

T = ^T b ~ c T M t m t- cd + ^T bTi T d ~ c t m t- cd , (214) 

where T with the upper-case indices is the contravariant inverse of one with the 
lower-case ones, L _1 £ Herm (W <8> W), i.e., analytically, 

r°*r a = s\, T m T b ~ c = &i. (215) 

The metrics underlying (|214p are essentially curved and imply in a strong, 
non-perturbative nonlinearity. The corresponding action functional will be de- 
noted by 

I[T] = I[T, m] + I[T, a] = J Tdt. (216) 
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Its terms controlled by constants A, B were here denoted respectively by I[T, m], 
I[r,a]; the labels m and a refer respectively to "main" and "additional". The 
reason is that evidently the £?-term is a merely auxiliary correction and the 
model with A = would be meaningless. The A-term is a proper dynamics. 

After some calculations one finds the following expressions for variational 
derivatives: 



r , = -ir r a -r s! r'T a r, (217) 



"(f ffl -i\ fc r fe T H )r fca . (218) 



In analogy to doubly invariant Hamiltonian systems on semisimple Lie groups, 
one can show that the corresponding Euler-Lagrange equations 

™ -0 (219) 



sr m (t) 

are solvable in terms of the matrix exponential function 

r Ts (t) =G Tz eMEt) Z s, (220) 

where E £ L(W) ~W® W* and G — T(0) G Herm (w* ® W*^j is the initial 

position in the configuration space Herm (w* (8 W*^j . It is clear that (|220p is 

a solution of (j2 19[) for any G and E. One can also easily show that T(t) persists 
to be Hermitian for all t € R if E is G-Hermitian, i.e., if the sesquilincar form 

gEts '■= Gt z E z s (221) 

is Hermitian. 

This procedure is equivalent to the following one, based on the multiplication 
by matrix exponents on the left: 

T ¥s (t) = exp(FtyGz S , (222) 

where again G is arbitrary and F £ L (w ~ W (E> W must be G-Hermitian 
if r is to be Hermitian for any feK. In other words, the sesquilinear form 

W* <g> W* 3 {F G )- s = F-~G- S (223) 

must be Hermitian. 

It is seen that, depending on the choice of E or F, T may be oscillatory, 
exponentially growing or exponentially attenuating. All these situations may 
have something to do with decoherence, reduction, and other "paradoxes" of 
quantum mechanics. 

Remark: We were dealing here with the pure dynamics for the "scalar prod- 
uct" r, without any interaction with the "wave function" %j>. By the way, 
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the dynamics for T was purely amorphous in the sense that no fixed metric 
G £ Hcrm (w* ® W*^j was assumed. In principle we might assume some and 

admit for T [T] something similar to (fTTTj) . (T72]) . ([150]) . No doubt, such expres- 
sions are definitely less convincing from the point of view of first principles of 
symmetry. On the other side, without using any fixed G, we have at our disposal 
only the above geodetic models (1214j) for the pure dynamics of T. When some G 
is distinguished, then just like in (|210p we can construct potentials invariantly 
built of some basic invariants of the following form: 

Tr( G rf), p=l,...,n, (224) 

where 

G T r s := G^IW (225) 

Incidentally, let us mention that instead of I\, /2-controlled terms in ()171[) . (I172p . 
(|180|) one can try to use some simpler expressions based on the fixed G, just 
quadratic in generalized velocities with constant coefficients, although showing 
weaker symmetries, e.g., 

^QbcQdaf,^ + E. G ba G dc t _ bt _ d (226) 

Proposition: After introducing the new dynamical term (|214[) (perhaps with 
some mentioned modifications), we must go back to all previous Lagrangians 
for ip and modify the corresponding variational derivatives on two levels: 

1. One of subsystems of equations of motion has the following form: 

51 \T,i/>] 



sr ¥s (t) 



0, (227) 



where / is built of the sum of all possible Lagrangians. We have just 
calculated (|2l"7l) , ([2T8| for / [r]. But we must have (|227| just for the total 
I [r, i/j] built of the total Lagrangian L [T, -0]. Fortunately, in all previous 
Lagrangians T enters in a purely algebraic way, so the calculations are 
relatively simple. 

2. For all previous terms one must revise the variational derivatives 

ffi^O, ^ = (228) 

sr si? 

taking, however, into account that T is a dynamical quantity and the 
operation d/dt in Euler-Lagrange expressions 

dL d dL dL d dL , nnn . 
(229) 

introduces some new terms involving dT/dt. 
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2.4 Total system of essentially nonlinear dynamical equa- 
tions 

Let us review the corresponding equations of motion [18]. They are strongly 
nonlinear in (tp, T) and this nonlinearity is essential and non-perturbative. We 
have yet neither rigorous nor qualitative solutions, however the above remarks 
concerning the pure dynamics for T seem to indicate that this kind of nonlin- 
earity may be an alternative description of the open quantum system for ip with 
the surrounding symbolically represented by the dynamical T. The nonlinear 
interaction with T might be perhaps a good candidate for explaining the afore- 
mentioned " paradoxes" . This might be perhaps an approach alternative to that 
developed by Ingarden, Jamiolkowski, Kossakowski [T51 Q31 [T5] and others. 
Proposition: Let us consider the most general Lagrangian: 

L = ai iT- ab - ^ty 6 ) + a 2 T- ab t a ^ h + [aiT-ab + a 5 H- ab ] ^ b 

+ a 3 [r ba + a 9 ?v b ] r-ab + n[i>, r] deba r ab f - cd - v (v, r) , (230) 

where 

fity>, T] deM = a 6 [T da + a 9 ^ a iP d ] [T b5 + a 9 ^ b ] + a^^ 5 ^ 

+ a 7 [T ba + a 9 ?V b ] [L de + a 9 #Vl = r] fcad5 , (231) 
and the potential V can be taken, for instance, in the following quartic form: 

V(tP,T) = K(r- ab ^J b ) 2 ■ (232) 

The first and second terms in (|230p (those with a\ and a 2 ) describe the free 
evolution of wave function ip while T is fixed. The Lagrangian for trivial part 
of the linear dynamics (those with ay) can be also taken in the more general 
form / (Tabip a ip b ) , where / : K — > M. The term with as corresponds to the 
Schrodinger dynamics while T is fixed and then 

H\ = T a£ H- cb (233) 

is the usual Hamilton operator. If we properly choose the constants a.\ and a$, 
then we obtain precisely the Schrodinger equation. The dynamics of the scalar 
product r is described by the terms linear and quadratic in the time derivative 
of r. In the above formulae %p a = ip a denotes the usual complex conjugation 
and ai, i = 1,9, and k are some constants. 

Remark: The connection of the new constants ai in (|230|) with the previous 
ones, i.e., a, /?, 7 (|35|) . A and B (|214p . is as follows: 

A B 

a x =CL, a 2 = /3, a 5 = -7, a 6 = — , a 7 = — . (234) 

Applying the variational procedure we obtain the equations of motion as 
follows: 

S L 

-=- = a 2 Tab'ii' b + (o^rab - 2aiir a 6 ) ip b - 2astab'^ b Vcd'ip c ip d 
dip a ^ ' 
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and 



SL 



ab 



2a g (a 6 r ad r c - fc + a 7 r a6 f g£i ) v b (r dE + a 9 ^ d ) 

{2nT- cd i>^ d - a 4 ) T- ab ~ a 5 H ab - [a 3 a 9 + an] f B J V> h = 0(235) 



2ft[i/;,r] bad T 5(i + 2f2[i/;,r] bade f erf + (2nT- cd ^ d - a A ) ^ b 

2T da [a 6 T bs (r^ + a 9 W c ^ f ) + a 7 T bE (T^ + a 9 W¥)] t sd t gf 
a 2 % a ^) b + [a 3 a 9 + M ^ a V* + [a 3 a 9 - a x i] 4> a ^ b = 0, (236) 



where 
£i[-0,r] bade 



+ a 6 a 9 \tp & i) d + ^ a ?H [T bE + a 9 ^ b ] 

+ a 6 a 9 \Pi> b + lp £ ip b ] [r da + a 9 W4> d ] 

+ a 7 a 9 \W a i> b + ^j/ i ip b ] [T ds + a Wl>' t \ 

+ a 7 a 9 \Pii d + W4> d \ [r ba + a 9 ?> b ] (237) 

- a 6 \T ds Tf a (r 6E + a 9 J/*p b ) + r b( T' E (T da + a 9 ^ a ^ d )] f g/ 



- a 7 [r be T /a (r dE + a^^) + r de T /E (r ba + a 9 ^ s ^)] r s/- 

The Legendre transformation leads us to the following canonical variables: 
dL 

n b = —r- = a 2 T- ab W + anT^, (238) 
dip" 

dL 

7f a = — = a 2 T- ab ^ b - c^lW, (239) 
BT _ 

^— = a 3 [T ba + a 9 r^ b ] + 2Sl[iP, T] bad5 t 5d . (240) 



r a b 



til> 



The energy of our n-level Hamiltonian system is as follows: 

— - 81 81 81 — 

E = r—+i; b -^+t ab -. L = -(a i T- ab + a 5 H- ab )^ b 

dipa dtp" dT ab 

+ a 2 r ab 1p a ip b + nty, r] abEd r ab t 5d + k (iw> b ) 2 . 

Inverting the expressions (|238p . (|239p . (|240[) we obtain that 



a 2 a 2 



a 2 



a 2 



a b 



(241) 

(242) 
(243) 



45 



where 



a 7 



2 [V,r] 



^6 («6 + na 7 



ab 



T- ad Y- cb ^ d , 



1 + agO^T] 
a e + (n-l)a 7 / Oi[^,T 



a 6 {a 6 + na 7 ) \ 1 + a$6\ [ip, T] 
and then the Hamiltonian has the following form: 



H 



1 



-1 „ab„cd a 3, 



«4 



Of 

a 2 



(244) 

(245) 
(246) 

(247) 
(248) 

v^v 6 



9 

/-y^ o 

+ -fsity, r]-^ [r ba + ^ V] [r de + a 9 v c > d ] + « (IWV 6 ) . (249) 

Remark: Let us notice that if we suppose that the scalar product T is fixed, 
i.e., the equations of motion are as follows: 

a 2 V> - 2ai^ a + (2k0i [V>, T] - a 4 ) ^ Q - a 5 H a b i; b = 0, (250) 



then taking all constants of the model to be equal to except of the following 
ones: 

a 5 = -1, (251) 



ai = 



2' 



we end up with the well-known usual Schrodinger equation: 

ih^ a = H a b i[> b . 



(252) 



Remark: The first-order modified version of the Schrodinger equation is ob- 
tained when we suppose that T is a dynamical variable and a 2 is equal to 0, 
i.e., 



ihip a = H a b ip b 



ih 

— + a 3 a g 



r°T c - 6 v 6 



+ {2k0! [%l>, V] - a 4 ) ^ - 2a 8 r a T a i/> fc f sd ^ d 

2a 9 T a5 (a 6 t 5d t gb + a 7 t- cb t gd ) ^ (T ds + a 9 ^ d }253) 
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20[^,r] 6sd T 5 



cd 



ih 

— - a 3 ag 



ih 

— + a 3 a 9 



- 2a 6 r da r bs (r' E + « 9 ^V) t- cd t sf 
_ 2aiT dd T b-c p fs + ag ^/) f - cA t gf 

- {2k6x ty, r] - 04) - 20[V, r] bad5 t 5d . 
We can rewrite (|253[) in the following form: 

ihi> a = H eS a b ip b , 

where the effective Hamilton operator is given as follows: 

ih 



(254) 
(255) 



tl c S b — tl b — 



a 3 a 9 



cb 



+ (2rt0i [if>, T] - a 4 ) S\ - 2a 8 r a T eb f - ed ^ d 

- 2a 9 r ae (a 6 t sd t eb + a 7 t 5b t sd ] (T ds + a 9 ^ d ) 



2.5 Invariance properties of our general Lagrangian 



(256) 



So, if we investigate the invariance of our general Lagrangian (|230p under the 
group GL(n,C) and consider some one-parameter group of transformations 

{cxp (At) : r G K} , AeL(n,C), (257) 

then the infinitesimal transformation rules for tp and T are as follows: 



ip a h-> L a b Tp b , T aE 1 ► L a b L 5 s r bs 
where 



_ic_ t —Id 



(258) 



L a b - (5% + e^ a b , L- lQ b «5 Q b -e^ b , e m 0. (259) 

So, leaving only the first-order terms with respect to e we obtain that the vari- 
ations of ip and r are as follows: 

5i, a = eA a b ^ b , = eA a £ ^ 5 , (260) 

Sr ac = e ( A a bT bc + Ac_ T a^ t Sr _ b = _ e ( T _ b -Jc_ + T- ad A d b ) , (261) 



then 



i ^sr + ^ b 

e \ d ^d diP b 



= ra b (a 2 4> a + ai nP^A b d ^ d 

+ T- ab (u 2 ^ b -an^A a cV (262) 



and 



1 dL Ar 

-— ; OL ab 

esr ab 



a 3 (s b f + a 9 r- af ^») + 2fi [i>,r} badc r- af r- cd 



a 3 (6 a e + a 9 T gb ^ b ) + 20 [V>, if** r gb f 



cd 



A f b 

A s s (263) 
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If we consider some fixed scalar product T and take the r -hermitian A's, then 
A\ = T a ~ c A 5b , T e = A&Tf, l f = 2, (264) 
and therefore the above expressions can be written together in the matrix form: 

J (A) = Tr (Vl) , (265) 

where the hermitian tensor V describing the system of conserved physical quan- 
tities is given by the following expression: 

V = a 2 (V^ITo 1 + T^T^ + (aii - a 3 a 9 ) ^>t rr -i _ 2a3 r ( r 1 
- (a^ + a^T^T^ -2(T^ 1 Tlu[iP,T]+lu[iP,T]TT q 1 ) , (266) 

where 

c^,r? 5 = n[ib,rf ade r 5d . (267) 

Similarly for the To-antihcrmitian A's, i.e., when A^ — — A, we obtain another 
hermitian tensor W as a conserved value: 

J (A) = Tr (iWA) , (268) 

where 

iw = a 2 (^rro 1 - Tq ^v-v 1 ) + - a 3 a 9 ) v^rro 1 

+ + a 3 a 9 ) r 1 r^v t + 2 (r„ ^ ty, t}-lu [ii>, r] it- 1 ) . (269) 
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